LU T P 05-26 
|hep-lat/0506004| 

revised 
July 2005 



Masses and Decay Constants of Pseudoscalar Mesons to Two Loops 
in Two-Flavor Partially Quenched Chiral Perturbation Theory 



Johan Bijnens and Timo A. Lahde 



Department of Theoretical Physics, Lund University, 
Solvegatan 14A, S 223 62 Lund, Sweden 



Abstract 



This paper presents a first study of the masses and decay constants of the 
charged, or flavor-off-diagonal, pseudoscalar mesons to two loops for two flavors 
of sea-quarks, in Partially Quenched Chiral Perturbation Theory (PQ^PT). 
Explicit analytical expressions up to 0(p 6 ) in the momentum expansion are 
given. The calculations have been performed within the supersymmetric for- 
mulation of PQxPT. A numerical analysis is done to indicate the size of the 
corrections. 



PACS: 12.38. Gc, 12.39.Fe, 11.30.Rd 



Masses and Decay Constants of Pseudoscalar Mesons to Two Loops 
in Two-Flavor Partially Quenched Chiral Perturbation Theory 



Johan Bijnens 1 and Timo A. Lahde 1 

department of Theoretical Physics, Lund University, 
Solvegatan 14A, S 223 62 Lund, Sweden 

This paper presents a first study of the masses and decay constants of the charged, or flavor- 
off-diagonal, pseudoscalar mesons to two loops for two flavors of sea-quarks, in Partially Quenched 
Chiral Perturbation Theory (PQxPT). Explicit analytical expressions up to 0(p ) in the momentum 
expansion are given. The calculations have been performed within the supersymmetric formulation 
of PQxPT. A numerical analysis is done to indicate the size of the corrections. 

PACS numbers: 12.38.Gc, 12.39.Fe, 11.30.Rd 



I. INTRODUCTION 

At the present time, the most promising way to derive 
information on low-energy hadronic observables directly 
from QCD is by means of numerical Lattice QCD simu- 
lations. However, it is well known that such simulations 
are computationally very demanding if dynamical sea- 
quark effects (unquenched Lattice QCD) are to be taken 
into account. Although much progress has been made re- 
cently, simulations with sea-quark masses that are close 
to the physical u, d quark masses of a few MeV are still 
not available, whereas sea-quark masses down to several 
tens of MeV are becoming possible. The physical values 
of the quark masses then have to be reached by extrap- 
olation from these simulations. 

A generalization of QCD where the valence quark 
masses are different from the sea-quark masses is pos- 
sible, and is called Partially Quenched QCD (PQQCD). 
QCD may then be obtained as a continuous limit from 
this theory, in contrast to the quenched calculations 
where no such smooth connection exists. The reason for 
doing PQQCD Lattice simulations is that the valence 
quark masses can be varied much more efficiently than 
the sea-quark masses, which allows for more extensive 
studies of the quark mass dependence of the various ob- 
servables. 

The generalization of xPT to sea-quark masses differ- 
ent from the valence ones, i.e. to the partially quenched 
case 0] , provides a practical method for the extrap- 
olation to physical quark masses from the lattice simu- 
lations of PQQCD. The quark mass dependence of par- 
tially quenched chiral perturbation theory (PQxPT) is 
explicit, and thus the limit where the sea-quark masses 
become equal to the valence quark masses can be taken. 
As a consequence, unquenched xPT is included in 
PQxPT, and the free parameters, or low-energy con- 
stants (LEC:s), ofyPT can be directly obtained from 
those of PQxPT 0. The need for PQxPT calcula- 
tions to next-to-next-to-leading order (NNLO) has also 
been demonstrated by Lattice QCD simulations 00. 

The calculation of the charged pseudoscalar meson 
masses and decay constants to one loop (NLO) in 
PQxPT has been carried out in Refs. 0,0,0, and first 



results for the mass of a pseudoscalar meson at two loops 
(NNLO) in PQxPT have already appeared, for three fla- 
vors of degenerate sea-quarks 0] . The decay constants of 
the pseudoscalar mesons are also known to NNLO 9] for 
three sea-quarks. A full PQxPT calculation of the pseu- 
doscalar meson masses for three flavors of sea-quarks is 
in progress [lOj. However, as many Lattice QCD sim- 
ulations are performed with two rather than three sea- 
quark flavors, such expressions for the pseudoscalar me- 
son masses and decay constants to NNLO are also called 
for. 

This paper presents the full calculation of the masses 
and decay constants of the charged, or flavor off-diagonal, 
pseudoscalar mesons in NNLO PQxPT, for two flavors of 
sea-quarks (n sca = 2). The results are characterized by 
the number of nondegenerate valence and sea-quarks, de- 
noted (i V ai and d S oa, respectively. For the decay constants 
of the charged pseudoscalar mesons, the maximum num- 
ber of nondegenerate valence quark masses is cZ va i = 2. 
The degree of quark mass degeneracy in each result is 
sometimes also referred to with the notation d va i + d selL . 
The decay constant of the charged pion in the SU(2) 
symmetric limit for both valence and sea-quarks thus cor- 
responds to the 1 + 1 case. 

The analytical expressions for the NNLO shifts of the 
meson masses and decay constants are in general very 
long, but the expressions simplify considerably when the 
sea or valence quark masses become degenerate. In view 
of this, the NNLO loop results are given separately for 
each case of d va i + d sca considered. In the next sections, 
a short technical overview of the NNLO calculations is 
given, with an emphasis on the parts not covered in the 
previous publications 00], along with the full results for 
the NNLO contributions to the meson masses and decay 
constants. Finally, a numerical analysis of the results is 
given, along with a concluding discussion. 



II. TECHNICAL OVERVIEW 

The technical aspects of the PQxPT calculations to 
NLO have been thoroughly covered in Refs. 0,0. The 
new parts needed for NNLO are treated in Refs. . 
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For the sake of completeness, some of the issues which are 
directly relevant for the present two-flavor calculations 
are repeated here. 

Most significantly, the Lagrangians of PQxPT at 
0(p 4 ) and 0(p 6 ) can be directly obtained from the cor- 
responding Lagrangians of unqucnchcd rip flavor xPT 
by replacement of traces with supertraces. This follows 
because the derivation of these Lagrangians in Ref. 
only used relations which remain valid when all traces 
are replaced by the corresponding supertraces. For the 
definition of the supertraces in PQxPT see Refs fl!!. 
We work here in the supersymmetric version of PQxPT 
without the singlet <£>o as discussed in Ref. Q ■ The diver- 
gences of this version of PQxPT also follow directly from 
those of nj flavor \PT calculated in Ref . [T^| . provided 
that n/ is set equal to the number of sea-quark flavors. 
This is so since all the manipulations in Ref. |12| for nt 
flavors remain valid when traces are replaced with super- 
traces. Alternatively, this derivation of the Lagrangians 
of PQxPT can be argued for with the Replica method of 
Ref. jBj. 

A. Quark Masses and Low-Energy Constants 

The number of independent low-energy constants 
(LEC:s) in unqucnchcd and partially quenched xPT is 
slightly different, but the former are always linear com- 
binations of the latter. The number of LEC:s is also 
dependent on the number of flavors. Table [i] gives an 
overview of the various sets of LEC:s and the number of 
parameters up to NNLO for the different versions of xPT 
and PQxPT. 

TABLE I: The relevant sets of LEC:s, where the i+j notation 
indicates the number of physically relevant (i) and contact (j) 
terms in the respective Lagrangians. At NLO, the latter ones 
are usually denoted hi for rif = 2 and iJJ for n/ = 3. The 
relations between the various sets of LEC:s is discussed in 
the text. For simplicity, the superscript (npq) of the partially 
quenched LEC:s has been suppressed in most of this paper. 



n f 


XPT 
2 


X PT 
3 


X PT 

n 


PQxPT 
2 


PQxPT 
3 


LO 


F,B 


F ,B 


p (n) p(n) 


F,B 


Fo,B 


NLO 

i+j 


H 
7 + 3 


LI 
10 + 2 


11 + 2 


11 + 2 


11 + 2 


NNLO 

i + j 


53 + 4 


ci 

90 + 4 


112 + 3 


112 + 3 


112 + 3 



At lowest order, the LEC:s of PQxPT with n fla- 
vors of sea-quarks are the same as those of n flavor un- 
quenched xPT- At NLO with n flavors of sea-quarks, 



they are L r ^ npq ^ through L^™ 2 * 3 '. The two possible ex- 
tra terms li r ^ lpq } and can be removed using field 
redefinitions [Til Il2|. For the case of three sea-quarks, 
the L\ are simple linear combinations of the jj^ pq ^ ) 
and these relations can be found in Refs. 0, Sim- 
ilarly at NNLO, the C\ are linear combinations of the 

K^ 3pq \ and the corresponding relations can be found 
in Ref. [ll|. For three flavors of sea-quarks, one simply 
takes L r ^ pq) = L r £ pq) = 0. 

The unquenched two-flavor Lagrangian as defined in 
the second paper of Ref. 3] differs by an L T X1 type term. 
At NLO, this makes no difference since that term does 
not contribute, but in order to get the correct values at 
NNLO for the Z[,c[ and the L r ^ vq \ Kl (2pq) one should 
take L[ {2pq) = -11/4 and L[ {2pq) = 0. The relations 
between the NLO LEC:s are given in Table ITT1 and simi- 
lar linear relations between the NNLO LEC:s K^ 2pq ' of 
PQxPT and the of unquenched two-flavor xPT may 
be found in Ref. [llj . The equations given in the remain- 
ing sections of this paper only deal with n/ = 2 PQxPT, 

and therefore the superscript (2pq) of the Jj^ 2pq) and the 
j^r(2pq) nave j-, een suppressed in most of the formulas in 
the following sections. 

TABLE II: The correspondence between the NLO low-energy 
constants L\ of two-flavor PQxPT and the l\ of standard two- 
flavor xPT- The relations relevant for the NNLO low-energy 
constants K\ are given in Ref. [TH . 



n f = 2 X PT 


n s = 2 PQxPT 


jr 
£ 1 

ir 
L 2 
ir 

h 

jr 
£4 

ir 
L 5 
ir 
L 6 

jr 
L 7 


-2L r i2pq) +4L r 1 {2pq) +2L r 3 {2pq) 
4Ll (2pq) +4L r 2 (2pq) 
-8 L r 4 (2pq) - 4 Ll} 2pq) + 16 L r 6 {2pq) + 8 L r s (2pq) 
SL\ (2pq) +ALl (2pq) 

r r(2pq) 

—2 lf ( - 2pqS> 
-W L r 7 (2pq) - 8 Ll (2pq) 



The input quark masses rrii enter into the calculation 
in terms of the lowest order squared meson masses \ii 
which arc defined as usual in xPT, by x« = 2i? m^. In the 
present calculations, there are two valence inputs XI7X2, 
two sea inputs X3iX+ an d two ghost inputs X5;X6- In 
order for the disconnected valence quark loops to be can- 
celed, the masses of the ghost quarks are always equal 
to those of the corresponding valence quarks, such that 
X5 = Xi and \6 = X2- Explicitly, for d val = 1, we have 
Xi = X2 and for d va \ — 2, we have xi 7^ Xi- Similarly, for 
the sea-quark sector d sea = 1 implies X3 = X+ while for 
d sca = 2 we have X3 7^ Xi- The most general results (in 
the quark masses) for n sca = 2 are therefore those with 
d va i = 2 and d sca = 2. 



3 



B. Chiral Logarithms 

The finite parts of the one-loop integrals which are en- 
countered in the NNLO calculation of the masses and 
decay constants of the pseudoscalar mesons may be ex- 
pressed in terms of the chiral logarithms A, B and C . Of 
these, the logarithms A are defined as 

A{Xi) = -^i6X 4 log(xi/M 2 ), (1) 
A( Xi ,e) = A( Xi ) 2 /(2Tr 16 Xi) 

+ 7r 16Xl (^7l2 + l/2), (2) 

where ir^ = l/(167r 2 ). Note that the chiral logarithms 
which depend on e are naturally generated by the process 
of dimensional regularization and can be reexpressed in 
terms of the standard logarithms A and B. However, 
they have been retained in the final results for the sake 
of notational efficiency. The finite parts of the logarithms 
B are 



B{Xi,Xj,0,k) 



-7r 16 (l + log( Xl //J 2 )), (3) 
Xi log(Xi/M 2 ) - Xj log(Xj/M 2 ) 

- Tl6 , 

X l - Xj 

A(x 1 )B(x 1 ,X 1 ,0)/(tt 16Xi ) (4) 

-A{ Xl f/{2n wX l) 

+ tt 16 (tt 2 /12 + 1/2), 

XiB(xi,Xj,0)+A( Xj ), (5) 



where the B(xi, Xj > 0) integrals, which are symmetric 
under the interchange of Xi and Xj j have been introduced 
in order to make the symmetries in the end results more 
explicit. Finally the remaining one-loop integral C is 
given by 



C(Xi»Xi)X»»0) = - 7I "i6/(2xi)- 



(6) 



In all of these cases, the dependence on the subtraction 
scale fj, has been moved into the loop integrals. 



C. Two-loop Integrals 

The finite two- loop integrals H F , H F , B F X that appear 
in the NNLO calculation of the meson masses and de- 
cay constants may be evaluated using the methods of 
Ref. Note that the corresponding primed integrals 

H F , H[ , H^i indicate differentiation with respect to p 2 
and appear only for the decay constants. The notation 
used for the H integrals in this paper is otherwise similar 
to that of Ref. [14|. except that an extra integer argu- 
ment now indicates the propagator structure, such that 
e.g. 

H F {Xt,Xj,Xk,p 2 ) -> H F (n, Xl , Xj ,Xk,p 2 )- (7) 

The case of n = 1 indicates that the integral consists of 
single propagators only, as in Ref. 01 1 whereas n = 2 



indicates that the first propagator appears squared and 
71 = 3 that the second propagator appears squared. The 
cases with two double propagators that can appear in 
the calculations are n — 5, for which the first and second 
propagators appear squared, and n = 7 for which the sec- 
ond and third propagators are squared. Explicit expres- 
sions for n = 1 can be found in Ref. 0], and the other 
cases may be obtained by differentiation with respect to 
the masses of those expressions. 



Notation for Meson Masses and Decay 
Constants 



D. 



As explained in great detail in Refs. U,l?j,|9|, the neu- 
tral meson propagators in PQxPT have a complicated 
structure. In order to obtain a result which can be ex- 
pressed in terms of the chiral logarithms and two-loop 
integrals outlined in the previous subsection, that prop- 
agator has to be treated by partial fractioning, which 
produces expressions that consist of sums of single and 
double poles. The residues, which are rational functions 
of the sea and valence quark masses, pile up and pro- 
duce very long and awkward end results. With a major 
effort, these can be simplified and in some cases the ex- 
pressions have been compressed by more than an order 
of magnitude. 

The residues R which are naturally generated in the 
calculations are, for d sea = 2, the single-pole residues 
Rj k and the double-pole residue Rf . The residues sim- 
plify when pairs of sea-quark masses become degener- 
ate, and hence the number of indices in the single-pole 
residue depends on the degree of degeneracy in the sea- 
quark masses. The residues relevant for the present two- 
flavor PQxPT calculations may be written in terms of 
the quantities R z a h , which are defined as 



R-ab — 

^abc 



Xa - 
Xa 



Xb, 
■ Xb 



R: 



Xa ~ Xc 

(Xa - Xb)(Xa 



Xc) 



bed 



Xa - Xd 



(8) 



and so on. Thus R z a b has the same dimension as Xi f° r 
an even number of indices and is dimensionless for an 
odd number of indices. It should be noted that the R z 
notation may also appear independently in the final re- 
sults. In such cases, the R z have been generated by sim- 
plification procedures. For d SC a = 2, the residues which 
originate in the partial fractioning of the neutral meson 
propagator are 



R 



jk 
Rf 

R'; 



i3ijk i 



(34tt • 



371 



R 



4tt 



Rl 



(9) 



where the neutral pion mass in the sea-quark sector is 
simply given by Xtt = X34 = (X3 + Xi) / 2 - Because of 
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this trivial relation, the above residue notation is highly 
redundant, a fact which has been exploited in the sim- 
plification of the end results. The quantity R v , which is 
defined as 



nil 
n ijk 



(10) 



appears throughout the final results and has therefore 
been given a separate name. For the case of d sell = 1, 
all residues associated with the sea-quark sector reduce 
to numbers. Some nontrivial residues may still appear if 
the valence quarks are nondegenerate, according to 



r>l rtz 

-ft,- — tin ■ 



(11) 



For convenience, Rf is also used for d va i = 1 in order to 
maintain a notation similar to the more general cases. 



E. Summation Conventions 

The expressions for the meson masses and decay con- 
stants are symmetric, within the sea and valence quark 
sectors, under the interchange of quark masses, and may 
thus be compactified by the introduction of summation 
conventions which exploit these symmetries. The sum- 
mations in the valence quark sector are similar to those in 
the three-flavor work of Ref . , whereas those of the sea- 
quark sector are considerably simpler in the two-flavor 
case. 

The sea-quark sector of the results involves the sum- 
mation index s, which may appear for the quark masses 
Xi, and among the indices of the residue functions R. A 
simple example is Xs — X3 + Xi- On the other hand, 
if the term in question consists of products of integrals, 
residues and quark masses, then the summation sign al- 
ways applies to the term as a whole, such that e.g. 

A( Xs )R 1 vs Xs = MxiJR^Xi- (12) 

i = 3,4 

Any contribution which is written in terms of the sum- 
mation index s explicitly fulfills the required symmetry 
properties in the sea-quark sector. 

The summation indices p and q, which refer to the 
valence quark sector, are only needed for the case of 
(fval = 2 and then always occur for the squared valence 
quark masses xi an d Xi- If the index q is present, there 
will always be an index p and the resulting sum is over 
the pairs (p, q) = (1,2) and (p, q) — (2,1). If only p is 
present, the sum is over the indices 1 and 2. An example 
is 

A(x P )R p q „Xp = A( X i)i?Lxi + [1-2]. (13) 

Any contribution written in terms of the (p, q) notation 
is thus symmetric under the interchange of the valence 
quark masses xi an d Xi- 



III. TWO-FLAVOR PSEUDOSCALAR MESON 
DECAY CONSTANTS 

The decay constants F a of the pseudoscalar mesons are 
obtained from the definition 



(OK(O)|0 a (p)) =iV2p»F a , 



(14) 



in terms of the axial current operator A£. The diagrams 
which contribute to that operator at 0(p 6 ), or NNLO, 
are shown in Fig. ^ Diagrams at 0(p 2 ) and 0(p i ) 
also contribute to Eq. (|14fl via the renormalization of 
the pseudoscalar meson wave function <fi a (p)- The decay 
constants are functions of the valence inputs \l > Xi an d 
the sea inputs Xi-iXtt which are defined in terms of the 
quark masses through Xi = 2_Bmi, and of the quantities 
Xij — (Xi + Xi)/2j which correspond to the lowest order 
charged meson masses. Other parameters include the de- 
cay constant in the chiral limit (F) , the quark condensate 
in the chiral limit, via (qq) — —BF 2 , and the LEC:s of 

C(p 4 ) and 0(p 6 ), i.e. the L 
for simplicity are denoted by L\ and K, 



r{2pq) and the K r S 2m \ which 





FIG. 1: Feynman diagrams at 0(p 6 ) or two- loop for the 
matrix element of the axial current operator. Filled circles 
denote vertices of the £2 Lagrangian, whereas open squares 
and shaded diamonds denote vertices of the £4 and Ce La- 
grangians, respectively. 



The decay constants of the pseudoscalar mesons are 
given in the form 



F 



phys 



F 



4 6)V8 + g 

F 2 



F 4 



0( P 8 



(15) 



where the 0(p ) and 0(p 6 ) contributions have been sep- 
arated. The contribution from the 0(p 6 ) counterterms 

has also been separated from the remaining part 



r(6) 
'ct 
r(6) 



^loops' which consists mostly of the contributions from the 



? (6) 



also 



chiral loops of 0{p ). It should be noted that £ loops 
receives a contribution from the chiral loops and coun- 
terterms of 0(p 4 ), which originate in the expansion of 
Eq. i|14|) to 0(p e ). The superscripts (v) and (s) indicate 
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the values of d va j and d se3j , respectively. As expected, the 
infinities in all expressions for the decay constant have 
canceled. The appearance of 'unphysical' B logarithms 
is, in part, due to the partial quenching, and to the fact 
that the results have been expressed in terms of the low- 
est order masses rather than the full physical ones. Thus 
not all of them correspond to the quenched chiral loga- 
rithms which arc ill-behaved in the chiral limit. 



A. Results for d va i = 1 

The NLO result for cZ va i = 1 is short, and will only be 
given for d soa = 2. The result for d sca = 1 can readily be 
derived from that expression. The combined NLO result 

I 



(loops and counterterms) for d sca = 2, is 

5^ 12 = 8^X34 + 4^X1 + l/2A(xi.), (16) 

and the contribution from the 0(p 6 ) counterterms to the 
NNLO expression for the decay constant is 

<5^ )12 = 8#[ 9 x? + 16i^ XiX34 + 8K r 2lX 2 s 

+ 32K r 22 xl 4 +8K r 23X l (17) 

At NNLO, the chiral loops form, by far, the largest con- 
tribution to the decay constant. As a straightforward 
derivation of the results for d sca = 1 from the d sca = 2 
case is tedious and complicated, the expressions for these 
two cases are given separately below. 



Cops = *i6 L r [-4xiX3 + 2x?-X^] - 2n 16 L[xi tt 16 L r 2 [ X \ + 3 X l] 

+ ir 16 L r 3 [-3xiX3 + 5/2 X f- 1/2 xl] + [-23/96 X iX3 + 1/64 X \ - 47/192 x §] - 2>1L\L\x\Xs 

- 32LJ 2 X § - 8 if xl + SLlfxi + 6A(xi)L r [xi+Rf\ - *A( Xl )L[ Xl - WA( Xl )L r 2Xl 
+ 6A( X i)L r 3 [xi+Rf] - 2A(xi)L r 5X i ~ 1/2 A( Xl )B( Xl3 , Xl3 ,0) Xl3 + 1/8 A( Xl ,e)ir mX3 

- A( Xi3 )tt 16 [l/6xi + 2/3x3] ~ 8A( Xl3 )L r oXl3 ~ 20 A( Xl3 ) L r sXl3 - 8A( X i 3 )^X3 + 4A( X i 3 )^Xi 
+ A(xi 3 ,e)7ri6X3 - 24A( X3 )L r lX3 ~ 6A( X3 )L r 2X3 + 12 A{ X3 )L\ X3 + 3/8 A( X3 , e) Xs 

+ 6B( Xl , Xu 0)L r R d lXl + 6B( X i,Xi,0)L r 3 Rixi - 2B( Xl , Xl ,0)L r 5 R d lXl 

- l/8B(xi,Xi,0, £ )^ 16 ^xi ~ 16 S(xi3,Xi3, 0)^X13X3 - 8 B( Xl3 , Xi3, 0) L r 5 X \ 3 

+ 32B(xi3,Xi3,0)LgXi3X3 + 16 J B(xi3,Xi3,0)L^Xi 3 - 1/8 # F (1, Xi, Xi3, Xi3, Xi) X3 

- 3/8ff F (l,xi3,Xi3,X3,Xi)X3 + l/8 J ff F (2,xi,Xi3,Xi3,Xi)^Xi + 5/12 H F ' (1, X i, Xi, Xi, Xi) X? 
+ ff F '(l,Xi,Xi3,Xi3,Xi) [1/8X1X3 - 1/2 Xi] + 3/8^ F '(l,Xi3,Xi3,X3,Xi)XiX3 

+ 1/2# f '(2,xi,Xi,Xi,Xi)#iXi + 3/8i/ F '(2,xi,Xi3,Xi3,Xi)^Xi 
+ l/4i/ F '(5, X i,Xi,Xi,Xi)(^) 2 Xi - 2^f'(3,Xi3,Xi,Xi3,Xi)^Xi 
+ 3/8 (l)Xi>Xi3,Xi3,Xi)Xi + 9/8 (l,X3,Xi3,Xi3,Xi)Xi 

- 3/8^ F '(2,xi,Xi3,Xi3,Xi)^Xi, (18) 

Cops - Tie LI [-4 x.Xi - 3 xl + 2xi + 2x3X4] - 2tt 16 ^x 2 - *i*L r 2 [3 X 2 +Xi] 

+ ir 16 L r 3 [-3x.Xi-2x'+5/2 X 2 + 3/2x3X4] + tt 2 6 [-23/96 x *Xi ~ 1/3 xl + 1/64 xl + 17/192 X3 X 4 ] 

- 32L^xiX34 - 32Lf X34 - 8 L r 5 2 x l + 8L\\ X \ + 6 A( Xn ) L r X „ - 8A( Xv ) L[ X „ 

- 2A( X „)L r 2X7r + 6A( Xn )L r 3 R^ X7T + AA( X „)L\ X ^ - 2A( Xn ) L\ R\ x X i 

- l/AA( x „)B( Xls , Xlsl 0)Rl sXls + l/8A( X7r ,e)n 16 R c lX „ + 6A( Xl )L r [R c lXl +R d ] - AA( Xl )L r lXl 

- 10A(xi)^Xi + 6A(Xi)L r 3 [R c lXl +R(\ - 2 A( X i) L\ R\ Xl - l/AA( Xl )B( Xu , Xu ,0)Rl wXu 

+ l/8A(xi,e)7T 16 i?^X34 ~ 16A(X34)iiX34 - 4 A(x 34 ) £2 X34 + 8 A(x 34 ) £4 X34 

+ l/4A(x34,e) ^16X34 - A(xi s )tti 6 [1/4 X34 + 1/6 xi,] - 4A( Xls )L r oXls - 10A( Xls )L r 3Xls 

- 4A( Xu )Ll Xu + 2A( Xls )L r 5Xl + l/4A(xi s ,£)^i6 [X34 + Xs] + 6 B( Xl , Xl , 0) L r R* X i 
+ 6B( Xl , Xl ,0)L r 3 Rf Xl ~ 2B( Xl , Xl ,0)L r 5 Rf Xl - 1/8 B(xi, Xi, 0, e) tt 16 R* Xl 

- 8 B( Xls , Xls ,0) L\x3A X \s - 4 B( Xls , Xls , 0)LgxL + 16 B(xu, Xu, 0) L r 6 X34Xis 

+ 8B(xi»,Xi».0)iSxL + l/16i/ F (l,X-,Xi s ,Xis,Xi)Ki s X- - H F (l, XllXls , Xls , Xl ) [l/8R\ vXu 

- 1/16 i^xi] - l/4^ F (l,Xi3,Xi4,X34,Xi)X34 + l/16ff F (2,xi,Xi s ,Xi s ,Xi)^iXi 
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+ i/iH F '(i lX ,,x,,xi,xi)(Rii) 2 xi + i/2H F ' (i, x *,xuXuxi)RliR c ixl 

- H F \l,Xir,Xu,Xis,Xi) [1/4 Rlxx\ + 1/16 KuXtXi] + H F ' (1, X i,Xi,XuXi) [l/6x? 
+ l/4(i?J) 2 X ?] + H F '(l,xi,Xis,Xu,Xi) [1/8 i?LxiXi,- 1/2 + 3/16 i^x?] 

+ l/4H F '(l,Xi 3 ,Xi4,X34,Xi)XiX34 + l/2i? F '(2, X i,X T ,Xi,Xi)i?? 1 i?fx? 
+ l/2ff F '(2, X i,Xi,Xi,Xi)^iXi + 3/16tf F '(2, X i,Xi s ,Xi s ,Xi)i?fx? 
+ l/4ff F '(5,xi,Xi,Xi,Xi)(«i) 2 Xi - 1/2 Hf (1,Xtt,Xis, Xu,Xi) RuRlvXl 

- #f (ijXis^i^xijxO-RL^LTrXi - H i {3,xis,xi,xis,xi) Rixl 

- 3/16fl^.'(l,x*,Xi-,Xi.,XiXi.Xi + fffi(l,Xi,Xi.,Xi«,Xi) [3/8 ii^xi- 3/16 iZ^Xi] 
+ 3/4^(1 ,X34,Xl3,Xl4, Xi)x? - 3/16ff 2 F /(2 

,Xi)Xis,Xis!Xi)-RiXi- (19) 



B. Results for d va i = 2 

In general, the expressions for ri va i = 2 are longer than 
those for ri va i = 1, as the number of independent mass 
combinations that can appear in the integrals is signifi- 
cantly larger. The most general two-flavor expression for 
the pseudoscalar meson decay constant is the result with 
d va i = 2 and ri sca = 2. However, in contrast with the 
results for ri va i = 1, the limit where the sea-quark masses 
become degenerate is nontrivial even at NLO. In the no- 
tation of Eq. I|15|) , the NLO contribution for ri va i = 2 and 

risea — 2 is 

6^ 22 = SLIxm + 4L$xi2 + MXv) [1/4 R% 

- 1/8 R°\ + l/lA{ Xps ) - 1/8 A{x*)Kn 

- l/8B( Xp ,Xp,0)R d p . (20) 



The NLO contribution for ri va i = 2 and ri sca = 1 may 
be obtained from Eq. I|20|l by carefully taking the limit 
Xa — ► X3- Effectively this means dropping the A( X ir) 
term and replacing A(x ps ) — » 2A(x p3 )- Furthermore the 
remaining residues reduce such that R° — > 1 and BP qv — > 
RP. The contribution from the 0(p 6 ) counterterms to 
the NNLO expression is 

<5^> 22 = 4K[ 9X l + 16^ xi2X34 + 8K r 2lX l 

+ 32 i^ 2X 2 4 + 8^3X1X2, (21) 

from which the corresponding expression for ri va i = 2 and 
dsca, = 1 may readily be inferred. The remaining NNLO 
contributions with ri va i = 2 are, for ri sea = 1 and ri sca = 2, 



Cops = T16 LI [-2X1X2- 4X12X3+ 4 X? 2 -Xi] - 2^16^X 12 - ^16^2 [x?2 + *xt] 

+ n w L r 3 [-3/2xiX2-3xi2X3 + 4x 2 2 -l/2x3] + ^6 [-17/192 X iX2 ~ 23/96 xi 2 X3 + 5/48 X 2 2 

- 47/192x1] - 32L^xi2X3 - 32 Lf xl ~ 8 X \ 2 + 8^ 2 X 2 2 + A( Xp )tt w [1/48 Xq + 1/48 X12 
+ 1/16 x 3 - 1/12 R> X12 - 1/8 RP q x 3 ] + A( Xp ) L r Q [ Xp + 4 RP q Xp + R*\ + A( Xp ) L r 3 [5/2 Xp + R p q X P 
+ 5/2 J$ + A( Xp )L\ [x3-2^x 3 ] - A{ Xp )Ll [1/2x12 + ^X12] + A( Xp ) 2 [9/128 - 1/32 i?^] 

- A( Xp )A( Xp3 ) [1/16+ 1/12 BP] - A( Xp )A( Xq3 ) [3/16+ 1/12 i$ + l/8A(x P )A(xn) 

+ A( Xp )B( Xp , Xp ,0) [3/8 Xp - 1/8 R p q Xp~ 1/32 i^ + l/64i$ - A( Xp )B( Xq , Xq , 0) [l/32J^ 

- 1/64 i?^] - l/4A(x P ) J B(x P 3,Xp3,0)xp3 - l/8A( Xp )B(xi,X2,0)R q p x P 

+ 1/8 A( Xp )C(x P ,Xp,X P ,0) Rf } X P - l/8A(x P ,e)7r 16 [x p -R p q Xa] + ^(XpsKie [1/12 x p3 - 1/4 x g3 

- 1/4 X3 ] - 4A(x p3 )L r Q x P 3 - 10A( Xp3 )L r 3 x P 3 ~ 4A(x p3 )^X3 + 2A(x p3 )^Xi2 - l/8A( Xp3 ) 2 
+ A(x p3 )B(x P ,x P ,0) [1/8Xp-5/8x p3 ] - 1/16 A( Xp3 )B( Xq , Xq , 0) < + l/6A(xp 3 )S(xi,X2,0)x3 

+ l/3A(x p3 )S(xi,X2,0,fc) + 1/2 i(x p3 ,e) ^16X3 - A(x0^(xa) [1/64- 1/16 R\R\] - 4A(xi 2 )^Xi2 

- 10l(xi 2 )^Xi2 + l/8A(xi 2 ) 2 - l/2A(xi 2 )S(xi,X2,0,fc) + 1/4 A( X ii,e) n w X12 

+ l/4A(xi 3 )^(X2 3 ) - 24l(x 3 )^X3 - 6A( X3 )^X3 + 12A( X3 )L r 4 x 3 + 3/16 A( X3 )B( Xp , Xp , 0) Xs 

- 3/8A(x 3 )B(xi,X2,0)x3 + 3/8 A(x3,e) ^16X3 + B(x P , Xp, 0) tt 16 [l/96i?^ X p + l/32i?^x« 
+ 1/16^X3] + B(Xp,Xp,0)LS^pXp + 5/2B(x P ,Xp,0)i^pXp + B( Xp ,X P ,0)Ll [2 Xp X3 
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- 4^XpX3 + 3^ X3 ] + B( Xp , Xp ,0)L r 5 [ X 2 p -2R p qXlX2 -l/2R d Xl2 ] + B( Xp , Xp ,0) L r 6 [A X 2 3 

- SR q p XpXs] + ^B( Xp , Xp ,0)L r 7 (R d ) 2 + B( Xp , Xp ,0)L r 8 [2 x \ - 4 B* Xp ] + B( Xp , Xp ,0) 2 [l/8 R p R d pXp 
+ 1/128 (R d ) 2 } - l/8B( Xp ,Xp,0)B(xi,X2,0)R q p R d pX p + 1/8 B( Xp , Xp , 0)C( Xp , X P , X P , 0) (R d p fx P 

- l/8B( Xp , Xp ,0 7 e)TT W R d Xp3 - 8 B( Xp3 ,X P 3, 0)L r 4 X P 3Xs - 4 B( Xp3 , Xp3} 0) L r 5 x 2 3 

+ l6B( Xp3 , Xp3 ,0)L r 6Xp3X3 + 8B( Xp3 , Xp3 ,0)L r sX 2 p3 + l/64B( X i,Xi,0)B(X2,X2,0)i2?i^ 

- 8B( X i,X2,0)L r 7 R d R d - 4B( Xl , X2 ,0)L r s R d R d + AC( Xp , Xp , Xp , 0) L\ R d XpX3 

+ 2C( Xp , Xp , Xp ,0)L r 5 R d x 2 p - 8C( Xp , Xp , Xp ,0)L r 6 R d XpX3 - 4C( Xp , Xp , Xp ,0) L r 8 R d x 2 p 

+ H F (l, Xp , Xp ,Xi2,Xi2) [1/8 Xp~ 5/64xi2 + 1/16 i^Xi 2 ] + # F (1, X P , Xis, X23, X12) [-1/16 

+ 1/16 x,- 1/8 R p qX3 ] + ^ F (l,Xi,Xi2,X2,Xi 2 ) [1/32x12-1/8^X12] 

- l/8iJ i? (l,Xl2,Xl2,Xl2,Xl2)Xl2 - 3/8iJ i? (l,Xl3,X23,X3,Xl2)X3 

+ if F (2,x P ,X P ,Xi2,Xi 2 ) [1/16^X12-1/32^X12] + !Z F (2,x P ,Xi2,X 9 ,Xi2) [l/16i^ X i2 

- 1/32^X12] + l/8ff F (2,x P ,Xi3,X23,Xi2)^X P 3 - l/64^(5,x P ,X P ,Xi2,Xi 2 )(^) 2 Xi2 

- l/32 J ff F (5,Xi,X2,Xi2,Xi2)^Xi2 - # F '(l,Xp,X P ,Xi2,Xi 2 ) [l/8 x P Xi2 + 15/64 X \ 2 - 3/16 RP q X \ 2 
+ 1/16 (R p ) 2 X 2 2 ] + # F '(l,X P ,Xi3,X23,Xi2) [l/8x 2 2 + 3/8^XpXi2 + l/8^x 9 Xi2-l/8i?«Xi2X 3 ] 

- J ff F '(l,Xi,Xi2,X2,Xi 2 ) [1/32 x? 2 - 3/8 i?^?x? 2 ] + 1/8 H F ' (1, X i2, X12, X12, X12) Xi 2 

+ 3/8ff F '(l,xi3,X23,X3,Xi 2 )xi2X3 - # F ' (2, X P , Xp, X12, X12) [1/16 R p q R d pX \ 2 - 5/32 ^ X ? 2 ] 

- J ff F '(2,Xp,Xi2,X< z ,Xi 2 ) [l/16i?^x 2 2 -l/32^X 2 2 ] - 1/8 H F ' (2, Xp , Xl3 , X23 , Xl2 ) R d Xp3Xl2 
+ 5/64ff F '(5,Xp,Xp,Xi2,Xi2)(^) 2 X 2 2 + l/32iJ F '(5,xi,X2,Xi2,Xi 2 )i?^X 2 2 

+ ^f'(l,Xp,Xp,Xi2,Xi 2 ) [5/4x 2 2 -l/2^x 2 2] - ^f'(l,Xp,Xi3,X23,Xi 2 )^X 2 2 

- i/f'(l,Xp3,X 9 3,Xp,Xi2)x 2 2 + ^f'(l,Xi2,Xi,X2,Xi2) [l/4 x? 2 - 1/2 R\R\ X? 2 ] 

- 1/4 ijf (3, Xi2,Xp,Xp,Xi2)i?«i?pX 2 2 + l/4i/f'(3,xi2,Xp,X?,Xi2)^i?pX 2 2 

- l/8^f'(7,Xi2,Xp,Xp,Xi2) (i?p) 2 X 2 2 - 3/8ff 2 F 1 '(l,Xp,Xp,Xi2,Xi2)Xi 2 
+ 3/8ff£'(l , Xp,Xl3,X23, Xi 2 )^X 2 2 + 3/8<(l 

, Xp3, Xg3, Xp, Xl2)i??X 2 2 

- 3/8ff 2 i'(l,Xp3,X 9 3,X 9 ,Xi2)^X 2 2 + ^2 F i'(1,Xi2,Xp,Xp,Xi2) [15/64 x 2 2 - 3/16 X ? 2 ] 

- ff 2 F 1 '(l,Xi2,Xi,X2,Xi2) [3/32 x? 2 - 3/8 R\R 2 lX 2 12 ] + 3/8<(l, Xi2,Xi2,Xi2,Xi2)x? 2 

+ 9/8iJ 2 F '(l,X3,Xl3,X23,Xl2)x 2 2 - 3/8i7 2 F 1 '(3,Xp3,Xp,X 9 3,Xl2)-RpX 2 2 

- ^2i'(3,Xi2,Xp,Xp,Xi 2 ) [3/16i?^x 2 2 -3/32i?^x 2 2 ] - #£'(3, X i2, Xp, X„ X12) [3/16 i?«i^ X 2 2 

- 3/32i?^ X 2 2 ] + 3/64 J ff 2 F '(7,xi2,Xp,Xp,Xi2)(i?p) 2 Xi 2 + 3/32 ff 2 F '(7, X i2, Xi, X2, X12) R d R d 2 Xi 2 , (22) 

Cops = *16 £0 [-4XxXi2"3 X 2 -2xiX2 +4 x? 2 + 2X3X4] - 2^ 16 ^ X 2 2 - tt 16 [3 X 2 +Xi 2 ] 

+ tt 16 Lg [-3 X.X12 - 2 x 2 - 3/2 xiX2 + 4 x? 2 + 3/2 X3 X 4 ] + n 2 16 [-23/96 x.Xi2 - 1/3 X 2 - 17/192 X1X2 
+ 5/48 x? 2 + 17/192x3X4] - 32i^ X i2X34 - 32 if X 2 3i - 8L r2 X \ 2 + 8L r2 lX \ 2 
+ A( Xp ) tt 16 [-1/12 RP qn X12 - 1/8 R p q7r X34 + 1/48 R c p Xq + 1/48 R c p Xl2 + 1/16 i£ X34 ] 
+ A( Xp ) L r [4 R p q7T xp + ^ Xp + Rt] + MXr) LI [R p qn Xp + 5/2 R c p Xp + 5/2 R d ] 

- A( Xp )Ll [2^X34-^X34] - A(xp)£S [^xi2 + 1/2^X12] + A( Xp ) 2 [1/16+ 1/32 (^J 2 

- 1/32^^ + 1/128 (i?p 2 ] - A(xp)A( X p S ) [1/24^ + 5/48^-7/96^] 

- A(xp)A(x 9S ) [1/24^ + 1/16^ + 1/32^] - A( Xp )A( Xv ) [1/32^^2-1/64^^] 

+ l/8A(xp)A(xi 2 ) + l/12 J 4(xp)A(x 34 )i?;p + A(xp)B( X p, x P , 0) [1/4 Xp - 1/8 R p qjr R pXp - 1/32 i^i^ 
+ 1/8 (R c p ) 2 Xp + 1/64 R c p R d ] + l/8A( Xp )B( Xp , x ^,0)R^ p R c p Rl 12Xp - A( Xp )B( Xq , Xq , 0) [1/32 i^i^ 



- 1/64^] - l/8A( Xp )B( Xps ,x P s,0)R p SJT X P s - 1/8 A( Xp )B( Xl , X2, 0) R q pn R c p x P 

+ 1/8 A{ Xp )C{ Xp , Xp , Xp ,0) R c p R d Xp + Mx P ,e)Ki6 [-1/8 Xp + l/8R p qnX34 + l/8R; pX34 ] 
+ A(xp S )7Ti 6 [l/24 Xps - 1/8 x 9S - 1/8x34] - 2A(x P s)L r oXps - 5 A( Xps ) L r 3 Xps - 2A( Xps ) L r iX34 
+ A( X ps)L r 5Xl2 - 1/32 A{ Xps ) 2 + A(xps)A(x„) [7/96 R^ - 5/48 i£ 8 - 1/32 B? qq + 5/48 R* qs - 1/24 
+ A( Xps )B( Xp , Xp ,0) [1/16 iCxp- 5/16 i?^ XPS ] - ^(X PS )^(XP,X W ,0) [l/12 i^i?^ Xp 

+ l/ei^i^xp-] - l/32A( Xps )B(x q ,X g ,0)R d g + l/12A( X p S )B(xi,X2,0)RLXs 
+ l/6A( Xps )B{xuX2,0,k)R q sn + l/8A{ Xps ,e)ir 16 [ X34 + Xs] - 1/16 A( Xp3 )A( Xp4 ) 
+ l/16A(x p3 )A( Xq 4) + A(x,)tt 16 [1/24 < 12Xt + 1/16 < 12 xi 2 ] + A( x .) L r [6 R\ 2 X , + K12 X,} 

- 8A(x*)L r lX « - 2%)ik, + A( Xir )LS [6^ 2 x, + 5/2Ki2X,] + A( Xjr )L; [4** + K12X*] 

- A(xx)^ [2^x12 + 1/2^12X12] + 1/128 A( Xff ) 2 (i?^ 12 ) 2 + l/12A(x*)A(X3i)Kn 

+ A( x ^)B{ Xp , Xp , 0) [-1/8 RP q7r R; p Xp + 3/16 ^ Xp - 1/8 (^) 2 Xp - 1/16 R d p + 1/64 i^i^ia] 

- l/8A(x 7r )B(x p ,x.,0)^ 7r ^ 2 ^ g7r Xp - l/8A( x ,)B( Xps , Xps ,0)i?; sXps 

- l/16A(x.)B(x.,X-,0)Ki2X34 - 1/8 A(x„)B(xi,X2, 0)^X^X34 

+ l/8A(x w )(7(x P ,X P ,X P ,0)^Xx P + l/8A( XT ,e)7r 16 [ Xot - i?? 2 Xx - K12 Xn] 

+ A( Xl )A( X2 ) [-1/32 R p q7r R q + 1/16 R^R^ + l/eAR^] - 4 A( Xl2 ) L\ Xl2 - 10 1( X i 2 ) L£ X12 

+ l/8A( Xl2 ) 2 - l/2A( Xl2 )B( Xl;X2 ,0,fc) + l/4A(xi 2 ,e)7ri 6 xi 2 - 16 A( X34 ) £J X34 

- 4A(x34)i 2 X34 + 8A(x 34 )i4X34 + 1/8 A(X34)S(X P , Xp , 0) X34 

+ A( X34 )B( Xp , XlT , 0) [1/12 i? 9 %i?: i2 Xp - 1/6 i? 9 %Ki 2 X 34 ] - 1/6 i( X34 )S( Xp , X7r , 0, k) R^ p Rl 12 

+ 1/8A(X34)B(X„ X „0)R1 12 X34 + ^(X34)S(X1,X2,0) [-1/4X34 + 1/6^X12-1/4^X34] 

- 1/6 A( X 34)B(Xl,X2, 0,^)^2 + 1/4 A(X34,£) ^16X34 + l/16A( Xls )A( X2s ) 

+ B( Xp , Xp ,0)tt 16 [1/96 R d pXp + l/32^ Xg + l/16^X34] + 5(x P , Xp, 0) #p Xp 

+ 5/2 B(xp, Xp, 0) L r 3 R d p Xp + B( Xp , Xp , 0) L\ [-4 R p qn XpX34 + 2 ^ XpX34 + 3 R d p X34 ] 

+ B(xp, Xp , 0) [-2 R p qn X i X2 + R p X 2 p ~ 1/2 R d X i 2 ] + B( Xp , X P , 0) L r a [8 XpX34 - 4 XpX34 

- 4 ^ X34] + 4 S(x P , Xp, 0) (i?^) 2 + B(x P , Xp, 0) L£ [4 i?^ XlX2 - 2 i£ X 2 + 2 (i?^) 2 ] 

- S(Xp, Xp, 0) 2 [1/8 Xp - 1/8 R c p R d x P - 1/128 (R d ) 2 ] 

+ l/8B( Xp , Xp ,0)B( Xp , X „,0)R z qnp R d RZ 12Xp - l/8B( Xp ,Xp,0)B(xuX2,0)R%R d p Xp 

+ l/8B( Xp , Xp ,0)C( Xp , Xp , Xp ,0)(R d ) 2 Xp - l/16B( Xp , Xp ,0,e)7r 16 [R d Xp + R d X34 ] 

+ 8B{ Xp , Xv ,Q)L r 7 R qpn R p R z ^ 3ip + 4B( Xp , X tt, 0) L r s R qp7r R d Rl 34p - 4 5( Xps , Xps , 0) L 4 XpsX34 

- 2 5( Xps , Xps ,0)L£ X 2 s + 8B( Xps , Xps ,0)L£ XpsX34 + 4 B( Xps , Xps , 0) L r s X 2 S 

+ 2B( X7r , X7r ,0)^Ki 2 X 2 + B( X7r , X7r ,0)^Ki 2 x' - 4B( X7r , X7r ,0)^iC 2X 2 

- 4B( Xw , x ^0)^( J R: 4 ) 2 Ki 2 - B( X *,X*,0)L r 8 [K12 xl + K12 xl] 

+ l/64B( X i,Xi,0)B(X2,X2,0)^ - 8B(xi,x 2 ,0)^^ J R 2 d - 4 B( Xl) X2 , 0) L r 8 R d R d 
, Xp , Xp , °)^4- R pXpX34 + 2C( Xp 

, Xp , Xp , 0)LlR d pX 2 p - 8C( Xp , Xp , Xp , 0) L 6 Rp XpX34 

- 4C( Xp , Xp;Xp ,0)^^ x 2 + ff F (l,Xp,X P ,Xi2,Xi2) [1/8 X P - 1/16 X12- 1/16 (R p q7r ) 2 Xi2 
+ 1/16 R p q7r R c p X12 - 1/64 (^) 2 X i 2 ] + ff F (l, Xp , Xls , X2s , Xl2 ) [-1/8^x^ + 1/16^X12 

- 3/32 xp + 1/32 i?L x 9 + 1/16 ^ Xp S ] - 1/8 H F (1, Xp3 , X 9 4, X34, X12) X34 

+ H F (1, x., Xp, X12, X12) [1/16 i?^Ki2 X12 ~ 1/32 R c p R v wl2 X12} - 1/64 H F (1, X „ X„ X12, X12) (Ki 2 ) 2 

+ l/32i/ F (l, XT , Xls , X2s , Xl2 ) [-RZ s Xp + R ps X q + RZ ps X* + RZi2Xs] 

+ H F (1, xi, X12, X2, X12) [1/16 R p q7r R c q X12 - 1/8 X12 - 1/32 R^ Xl2 ] 

- l/8i/ F (l, Xl2 , Xl2 , Xl2 , Xl2 ) Xl2 + H F (2, Xp , Xp , Xl2 , Xl2 ) [1/16^^X12-1/32 7^X12] 

- l/32H F (2,Xp,X-,Xi2,Xi 2 )^Xi2Xi2 + i/ F (2,Xp,Xi2,X 9 ,Xi2) [1/16 R q pn R d X12 - 1/32 R c q R d X i 2 ] 
+ 1/16 H F {2, Xp , Xls , X2s , Xl2 )R d Xps - l/64tf F (5,Xp,Xp,Xi2,Xi2)0R p ) 2 Xi2 
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- l/32i/ F (5,xi,X2,Xi2,Xi 2 )^i?2Xi2 + H F '(l, Xp ,X P ,Xi2,Xi2) [-1/8 X pXi2 - 3/16 X ? 2 

- 1/16 {R^fx\ 2 + 3/16 R p q7r R c pX 2 12 ~ 3/64 {R c p ) 2 X \ 2 } + H F ' (l lXp , XlsiX2s , Xl2 ) [-5/16^^X12 
+ 1/8 RP qn Xp sXi2 + 3/32 RP„ Xq Xi2 + 7/32 B? n Xl2Xs - 1/16 R c p XpsXu] 

+ l/8i/ F '(l,X P 3,X 9 4,X34,Xl 2 )Xl2X34 + H F ' (1, X ,r, Xp, Xl2, Xl 2 ) [V 4 -Rq,r#12 Xl2 _ V 16 R qn R ll2 Xl2 

+ l/8^X?2 + V32^ U Xi2] + ff F '(l,X7r,X.,Xi 2 ,Xi 2 ) [l/4(i?? 2 ) 2 X? 2 + l/4i?r 2 K 12 X? 2 
+ 5/64« 12 ) 2 x y + # F '(l,X^Xi s ,X2 S ,Xi2) [-3/32 i?; sXp Xi 2 - 5/32 i?; sXg Xi 2 + 1/4 K 2 X? 2 

- l/32^ psX7 rXi2-l/32J^ 12 xi 2 Xs] + # F '(l,Xi,Xi 2 ,X 2 ,Xi 2 ) [-1/16 RP qn R q X? 2 + 3/8 Rl n R 2 ln x? 2 
+ l/32E^Xi2] + l/8i? F '(l,Xi 2 ,Xi 2 ,Xi 2 ,Xi 2 )x? 2 + ^ F '(2,X p ,Xp,Xi 2 ,Xi 2 ) [-l/16^ 7r ^x?2 

+ 5/32 R c p R d pX 2 12 ] + i? F '(2,x P ,X.,Xi 2 ,Xi 2 ) [l/8 X?2 + 1/32 J?Xi2 Xia] 

- H F '(2, Xp ,Xi 2 ,X 9 ,Xi 2 ) [l/16i?^ X ? 2 -l/32i?^X? 2 ] - 1/16 H F ' (2, Xp , Xu, X2s, X12) R d p X P sXi2 
+ 5/64^ F '(5,x P ,Xp,Xi2,Xi2)(^) 2 Xi2 + l/32ff F '(5, X i,X 2 ,Xi 2 ,Xi 2 )i?^ 2 X? 2 

+ H[' (1, x P , Xp, X12, X12) [X12 + 1/2 W 2 X? 2 - 1/2 R^Rp X? 2 + 1/4 (^) 2 X? 2 ] 
+ l/2i/f'(l,Xp,Xi s ,X 2s ,Xi 2 )^ 7r ^ g pX 2 2 - l/2i? 1 F '(l,Xp S ,X gs ,Xp,Xi 2 )i?Lx 2 2 

- l/2i/ F '(l,Xp S ,X gS ,X 7 r,Xl 2 )^ 2 ^sp ; rX 2 2 + l/4fff (l,Xl2,Xp,X7r,Xl2) [ RP qnKl2 Xu 

X12] - h[ (i,xi 2 ,X7r,X7r,xi 2 ) [i/4^ 2 -R^i 2 x? 2 + 1/8 (Ki 2 ) 2 x 2 2 ] 

,Xi 2 ,Xi,X2,Xi 2 ) [l/4 Rq^Rq X12 - l/ 2 ^^ X12] + l/4-H"f (3,Xi 2 ,Xp,Xp,Xi 2 ) [R P q7! R p Xi2 

- fi^X? 2 ] + l/4i?r(3,xi 2 ,Xp,X 9 ,Xi 2 )^i?pX 2 2 + l/4if F '(3,xi2,Xp,X^,Xi2) J Rr 2 ^ g ,r^pX 2 2 

- 1/8 H[ (7, Xi 2 ,Xp,Xp,Xi 2 ) (R p ) 2 Xi2 - 3/8 77^ (l,Xp,Xp,Xi 2 ,Xi 2 )Xi 2 

- 3/16 #£(1 X12 + 3/16^ (1 

, Xps, Xqsi Xp, Xl2) [Rqir Xl2 + #?tt Xl 2 - -Rp X12] 
+ 3/16 ff 2 F '(l, Xps, Xqs, Xq, Xl2) R^Rspq X 2 12 + 3/16 (1, Xps, Xgs, Xtt, X^) R^lq^Kjm X\2 

- 3/16i/ F '(l , X7r,Xis,X 2 s,Xi 2 )^r 2 ^si 7r -Rs 27r Xi 2 + #21 (l,Xi 2 ,Xp,Xp,Xi 2 ) [3/16xi2 + 3/16 (i?q W ) 2 Xi 2 

- 3/16 fl^x?2 + 3/64 {R c p f X \ 2 } + Hg(l,xi2,Xp,X*,Xi2) [-3/16 i?^Ki 2 X? 2 + 3/32 R c p Rl l2 x \ 2 \ 
+ 3/64 ff F (1 

, Xl2, Xtt, Xtt, Xl2 )(Ki 2 ) 2 x 2 2 +<(i ,Xi 2 ,Xi,X 2 ,Xi 2 ) [-3/16i?^i?gXi2 + 3/8 i?2 W i? 2 7r Xi 2 
+ 3/32i?^x? 2 ] + 3/8i/ F '(l,xi 2 ,Xi 2 ,Xi 2 ,Xi 2 )x 2 2 + 3/8 ff F '(l, X 34, Xp3, X 9 4, Xi 2 ) X? 2 

- 3/16i/ 2 F '(3,Xp S ,Xp,X 9S ,Xi 2 )#pX 2 2 - # 2 F i'(3,Xi 2 ,Xp,Xp,Xi 2 ) [3/16 R p qn R d x\ 2 ~ 3/32 i^x? 2 ] 

- <(3 , Xi2,Xp,X<?,Xi 2 ) [3/16 R p7I R p Xi 2 - 3/32 R c q R pX \ 2 ] + 3/32 H 21 {i, Xl2lXp , XlT , Xl2 ) R p R v nl2X \ 2 

+ 3/64i/ F '(7,Xi 2 ,Xp,Xp,Xi 2 )(i?p) 2 X 2 2 + 3/32i/ F '(7,xi 2 ,Xi,X 2 ,Xi 2 )^ 2 X 2 2 . (23) 



IV. TWO-FLAVOR PSEUDOSCALAR MESON 
MASSES 

The pseudoscalar meson squared masses Mp S are cal- 
culated from the resummation of the irreducible self- 
energies, according to 



M| s = Ml 



S(Mj 



ps)' 



(24) 



where Mq is the lowest order meson mass, and £ is the 
sum of the irreducible self-energy diagrams. The self- 
energy contributions — at 0(p e ), or NNLO, are shown 
in Fig. |21 The pseudoscalar meson mass is a function 
of the same mass inputs and LEC:s as outlined for the 
decay constants in the previous section, although the set 



of 0(p e ) LEC:s that contribute to the meson mass is 
slightly different from that of the decay constant, in par- 
ticular the number of 0(p 6 ) LEC:s is larger for the meson 
mass. The results for the pseudoscalar meson masses are 
given in the form 



PS 



Mi 



j(4)v< 



? (6> 



r(6)vs 
loops 



F A 



+ 0(p s ), (25) 



where the 0(p 4 ) and 0(p 6 ) contributions have been sep- 
arated. It should be noted that the lowest order mass 
has not been factored out of the results, as this is not 
meaningful for the results with e? va i — 2. Similarly to the 
treatment of the decay constant, the contribution from 
the 0(p 6 ) counterterms has also been separated from the 
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rest of the NNLO result, which consists mostly of the 
contributions from the chiral loops of 0(p 6 ), although 
there is also a contribution from the diagrams at 0(p 4 ) 
through the expansion of Eq. 1251) to 0(p 6 ). The super- 
scripts (v) and (s) indicate the values of <i va i and <i sea , 
respectively. As expected, the infinities in all expressions 
for the meson mass have canceled. 




FIG. 2: Feynman diagrams at 0(p r ') or two- loop for the irre- 
ducible self-energy contribution — iS to the pseudoscalar me- 
son mass. Filled circles denote vertices of the £2 Lagrangian, 
whereas open squares and shaded diamonds denote vertices 
of the £4 and C§ Lagrangians, respectively. 



dvai = 1 and d sca = 2 is 



£(4)12 = _ 16LlxiX3i -SLlxi 
+ 32LgxiX34 + lQL r 8 xi 

- 1/2 A( X ,)R^ Xi - l/2A{ X i)R c 1 Xi 

- 1/2 B( X i,Xi, 0)fl?xi, (26) 



from which the d sca = 1 case can be obtained by carefully 
taking the limit Xi ~ * X3- I n that limit the A(Xir) term 
vanishes, and the remaining residues reduce such that 
Ri — > 1, yielding the correct result for d va \ = 1 and 
dsca = 1- The tree level contribution from the 0(p 6 ) 
Lagrangian is 



<$ )12 = - 32K r 17 xl - 64^x1x34 ~ 16^ 9 x? 

- 32K r 2oX lx34 - l6K r 2lX iX 2 s 

- MK^xixli ~ IKK^xl + 48K r 25 xl 

+ K26 [l6xiX^ + 64 Xl X3 4 ] + 192K r 27 xixli 
+ 32^x1+64^x1x3^ (27) 



A. Results for d va i = 1 

The lowest order mass, M§ , is for d va i = 1 equal to xi- 
The NLO, or 0(p i ) contribution to the meson mass for 



from which the result for d soa = 1 may again be obtained 
by setting X 4 ~~ > X3- The remaining contributions to the 
meson mass at 0(p e ) we again give separately for d aca = 1 
and d sca = 2. These contributions are 



Sloops = me L r [2 X iX3 + 8XiX3-4 Xl ] + ^ w L r lX \ +n 16 L r 2 [6 X iX3 + 2xi] +n 16 L r 3 [ X i X l + 6 X \x3 - 5 xl] 
+ n 2 ie [47/96x1X3 + 35/144 X ? X3 - 17/288 X ?] + 256L^X?Xs - 512L^xiX3 - 256L^ X ?X3 
+ \2%L\L\ XX \X3 + 256Lf xixl - 256L^x?X3 - 128L^x! + 64LJL^ X f + 64 if X \ 

- 2hQLlL\ lX \x3 - 128LKa? - l2A( X i)L r [xl + Rfxi] + 8A( X i)L r lX i + 20A( X i)L r 2 xl 

- 12A( X i)L r 3 [xl + Rixi] + 16A(xi)^ X iX3 + A{ X i) L r B [16 X? + 8 Rf Xl] 

+ WA( X i)L r 6 [xi + Rfxi] + 32A( X i)L r 7 Rixi ~ 32l(xi)Lgxi + 5/8A( X i) 2 Xi 

+ A(xi)5(xi,Xi,0) [3/2x? + 1/4 Rixi] + 1/2 A( X i)C( X i, Xi, Xi, 0) Rf x? - A(xi, £ )7r 16 [7/4 x? 

- 1/4^X1] + 4/3l(xi 3 )7ri6XiX3 + 16 1(xis) £0X1X13 + 40A(xi 3 )£SxiXi3 - 32 A( X i 3 ) £5 X iXi3 
+ 64A(xi 3 )^8XiXi3 - -4(xi3) 2 Xi ~ 2 A(xi3)-B(xi,Xi,0)xiX3 - 2 A(xi3, e) Trie X1X3 

+ 48l(x 3 )£iXiX3 + 12 A( X 3) £2X1X3 - 48 A( X3 ) £4X1X3 + 48 A{ X3 ) L r 6 X iX3 
+ 3/4A(x 3 )£(xi,Xi,0)xiX3 - 3/4 A(x3,e) ^16X1X3 - 12 B( X i, Xi, 0) ££ Rfxl 

- 12B( X i,Xi,0)L r 3 Rfx 2 i - B( X uXi,0)L r 4 [24 X iX3 " 32 x?X 3 ] + 5(Xi, Xi, 0) ££ [4 X f + 24 R* xl] 
+ B(xi,Xi,0)L r 6 [32x1X3 - 48 X1X3] + 16 B( X i, Xi, 0) £' 7 (Ri) 2 Xl - B(xi,Xi,0) L r 8 [8xi+48i?fxi 
-8(^) 2 xi] +£(xi,Xi,0) 2 [l/2i?fx? + l/8(^) 2 Xi] + l/2B(xi,Xi,0)(7(xi,Xi,Xi,0)(^) 2 X 2 

+ 1/4 S(xi,Xi,0,e) ^e^xi + 16 C(xi,Xi,Xi,0) L^XiXs + 8 C(xi, Xl, Xl, 0) L r b Rf X \ 
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- 32<7( X i,Xi,Xi,0)L^x?X3 - 16 C( Xl , Xu Xu 0) L r 6 R d X \ + 5/6ff F (l, X i,Xi,Xi,Xi)Xi 
+ S" F (l,Xi,Xi3,Xi3,Xi) [1/4X1X3 -Xi] + 3/4H F (l,xi3,Xi3,X3,Xi)XiX3 

+ H F (2,xi,Xi,Xi,Xi)^Xi + 3/4H F (2, Xl , Xl3 ,xi3,Xi)i?fx? + 1/2 H F (5, X i,Xi,Xi,Xi) (Rffxl 

- 4ilf(3,xi3,Xi.Xi3,Xi)fliXi + 3/4^(1, xi,Xi3,Xi3,Xi)x? + 9/4^21 (l.X3,Xi3,Xi3,Xi)Xi 

- 3/4tf F (2, Xl , Xl3 , X i3, X i)i?f X 2 , (28) 



Cops = *i6 £5 [8x^X? + 6x'xi-4xiX3X4"4x?] + 4tt 16 ^ X ? + tti 6 L£ [6 X 2 X i + 2 X ?] 

+ 7r 16 L?J [6 X ^X?+4x'xi-3xiX3X4-5 X ?] + ^ 2 6 [35/144 XtX 2 + 2/3 x'xi - 17/96 X i X 3X4 

- 17/288 xl] + 256LJL5xiX34 - 512L^ X iX34 - 256L^x?X34 + 128 LJ^ X ?X34 

+ 256^ 2 xiX34 - 256L5LgxiX34 - 128 L^ X ? + 64L^ lX ? + 64L£ 2 X ? - 256L^ lX 2 X34 

- 128L^ lX ? - 12A( x .)iSKiX.Xi + 16A(x*)L r lX *Xi + 4 A( X „) L r 2 XvXl - 12 A( X „) L r 3 R^ X «Xi 

- 16A( X7T )L r 4 R c lXnXl + 8A( Xn )L r 5 [R^ X „ Xl + X \] + 16 A( Xn ) L r 6 B° X „ Xl 

+ 32A( Xn )L r 7 R* w34lXl - 32A(x*)LSi2r 1 Xi + 1/8 A( Xw f (R* u ) 2 X i + l/AA( Xw )A( Xl )R^R c lXl 

- l/3A( x .)i(x34)i?iiXi ~ 2/3A( Xn )A( Xls )R^ s R z UnXl - l/4A( Xn )B( Xv , Xv ,0)R? lXlXM 

+ 1/2 A(x*)B(xi,X«,0) RnKnXi + A( X% )B( Xl , Xl ,0) [3/4 R\ x \- 1/2 W) 2 X 2 - 1/4 Jifi2? X i] 
+ 1/2^( Xt )C(xi,Xi,Xi,0)^ 1 ^X? - l/4A(x.,e)7T 16 i?J X .Xi - ^A( Xl )L r [R° X \ + R d X i] 
+ 8A( Xl )L r lX l + 20A( Xl )L r 2X l - 12 A( Xl )L r 3 [R c lX \ + R d lXl ] + 16 A( Xl ) L\ R\ XlX34 
+ A( Xl )L r 5 [l6R c lX 2 + 8R d lXl ] + A( Xl )Lg [32 x 2 -16i^ XlX34 ] + 32 A( Xl ) U 7 R d X i 

- 32A( Xl )L^ x 2 + ^(xi) 2 [l/2 Xl + l/8(i?J) 2 X i] + 1/3 A( Xl )A( X3i ) R^ lXl 

+ 2/3A( Xl )A( Xls )Rl s Rl S7lXl + 1/2 A( Xl )B( Xl , Xw ,0) R^R c lX l + A( Xl )B( XuXu 0) [ x \ + 1/2 (Rl) 2 X 2 
+ l/AR\R d Xl ] + 1/2 A( Xl )C( Xl , Xl , Xl ,0) RIRfid - A{ X i,e)nw [l/2 X i X34 + 3/2 x 2 - l/4i?J XlX34 ] 

- l/2A( Xl3 )A( Xl4 ) Xl + 32A( X3i )L r lXlX34 + 8A( X34 )L r 2XlX3i - 32 A{ X3i ) L\ X i X34 

+ 32A(x 34 )i6XiX34 + l/2A( XM )B(x^X,^)RliX,Xi ~ A{ X3i )B{ Xl , Xv ,Q) [2/3^iX^Xi 

+ l/3i?? lX 2 ] + l/2A(x 34 )S(xi,Xl,0)^ 7r XlX34 - 1/2 A( X3 4,£) ^16X1X34 + A( X u) 7T16 [1/2 X 1 X 34 

+ l/6xiX-] + 8A( X i s )LSXiXia + 20i( Xls )^ XlXls - 16 A( X u) L\ XlXls + 32 A( Xu ) L r 8 XlXu 

- l/4A( Xls ) 2 Xl - A(xu)B(xi,x*,0) [2/3i?LxiXi s + l/3i?Lx 2 ] - A( Xu )B( Xl , Xl , 0) R] v XlXs 

- l/2A( Xls ,e)n w [ X i X3i + X i Xs ] + 8 B( X «, Xtt, 0) £^ i?^ X 2 X i + 4 £?( X7r , X7r , 0) R^ X \ X \ 

- 16B( X7r , X ^0)L^ lX 2 Xl - 16B( X7r}X7:7 0)L r 7 R^(RU) 2 X i - 4B( Xv , Xv ,0) L r 8 XlX 2 
+ R11 Xixl] + 32B{ XuX „,$)L r 7 R z r3A1 R d Xl + 16 B( Xl , Xw , 0) L r s R* 341 Rf X i 

+ 1/2 B(xi,X*,0)S(xi,Xi,0) R^Rf x l - 12B(xi,Xi,0)L5i?fx 2 - 12 B{ XllXll Q) L r 3 R d lX \ 
+ B( Xl , Xl ,0)L r 4 [8R c lX 2 lX34 + 24R d lXlX34 ] + B( Xl , Xu 0) L\ [4 R\ X \ + 24 R d x \] 

- 16B(xi,Xi,0)LS [i?J X 2 X 34 + 2^ XlX34 ] + 16B( X i,Xi,0)^(^) 2 Xi 

+ B(Xi,Xi,0)L5 [-8i?J X ?-48i?^ X 2 +8(i?^) 2 Xl ] + B( Xl , Xl ,0) 2 [1/2 x 2 + 1/8 (^) 2 Xl ] 

+ l/2S(xi,Xi,0)C(xi,Xi,Xi,0)(^) 2 X 2 + l/4S(xi,Xi,0,e)7T 16 i?f X 2 + 16 C( Xu X i, Xi, 0) L\ R d X \ X3i 
+ 8(7( X i,Xi,Xi,0)i^fx? - 32C( Xl , Xl , Xl ,0)Lgi?^ x 2 X34 - 16 C( X i, Xu Xu 0) L r 8 R d X \ 
+ l/2F F (l,x w ,X w ,Xi,Xi)(^u) 2 X 2 + ff F (l,X^,Xi,Xi,Xi)«n^X 2 

- if F (l, X „ Xls , Xls , Xl ) [1/2 Rl lX \ + l/8Rl lsXlxXl ] + ff F (l,Xi,Xi,Xi,Xi) [1/3 X 2 + 1/2(^) 2 X 2 ] 

- ff F (l,Xi,Xi„Xia,Xi) [3/4i?^ X 2 -3/8^ X 2 + l/8i?^ Xl ] + l/2F F (l,Xi3,Xi4,X34,Xi)XiX34 
+ F F (2,xi,X 7r ,Xi,Xi)^H^X 2 + ff F (2,Xi,Xi,Xi,Xi)^fx 2 + 3/8 tf F (2, X i, Xia, Xu, Xi) Rf Xi 
+ l/2F F (5,Xi,Xi,Xi,Xi)(^) 2 X 2 - #f(l, X ., X i„ X i s , X i)i?L^Xi 

- 2 J fff(l, Xls , Xls , Xl , Xl )^ s -RL 7 rX 2 - 2i/f(3, Xls , Xl , Xls , Xl )^ X 2 

- 3/8ff F (l, XT , Xls , Xls , Xl )K lsX 2 + Hf 1 (l,Xi,Xi.,Xi«,Xi) [3/4i?Lx 2 -3/8^ X ?] 

+ 3/2iJ F (l jX34,X13,X14, Xl ) X \ - 3/8i/ F (2, Xl , Xls , Xls , Xl )i?f x 2 . (29) 
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B. Results for d va i = 2 

The most general expression for the pseudoscalar me- 
son mass at n sea = 2 is the result which has d va i = 2 and 
4=ca = 2. In the notation of Eq. (|2*5|l . the lowest order 
mass Mq is equal to X i 2 . The contribution at NLO for 
rf va i = 2 and d sca = 2 is 

S^ 22 = - 16L r 4Xl2X34 - 8L r 5 x 2 12 + 32igxi2X34 
+ 16^xi 2 - l/2A( Xp )R p qn Xi2 
-1/2A( X „)RJ 2 X12, (30) 

from which the result for d va i = 2 and d sea = 1 is eas- 
ily obtained by taking the limit xa — ► X 3 such that 
RP„ -> i?£ and dropping the A(xn) term. The NNLO 



contribution from the 0(p e ) Lagrangian is given by 

5£ )22 = - 32*T[ 7X ? 2 - 64^[ 8X ? 2 X34 - 8K[ gX 2 Xl 2 

- 32^x^X34 - l&K r 2lX X2xl 

- 64 ^2X12X34 - 16^23X1X12X2 
+ 24^ 2 r 5 XpXi2 + [16*12x2 

+ 64 X ? 2 X3 4 ] + 192J^ 7Xl2X 2 4 + 32^ 9 x? 2 

+ 64^X12X34. (31) 

The result for c? va i = 2 and d sca = 1 may readily be 
inferred from Eq. (|31() by taking the limit xa —* X3- The 
loop contributions at NNLO, separately for d sca = 1 and 
rfsoa = 2, are 



Cops = ^16 ^5 [4xiXi2X2 + 2xi2X3 + 8x? 2 X3-8xy +47r 16 L[Xi2 + ^16 £2 [6 X12X3 + 2 X12] 

+ tt 16 L r 3 [3 X1X12X2 + X12X3 + 6 X? 2 X3 - 8 x? 2 ] + 7T? 6 [17/96 X1X12X2 + 47/96 X12X3 + 35/144 X ? 2 X3 

- 17/72 X12] + 256L^x? 2 X3 - 512i^ Xl2X ^ - 256i^ X ? 2X3 + 128L^ lX ? 2X3 

+ 256 if X i2X% - 2h6L\L\ X \ 2 Xz - 128L^5x?a + 64L^ lX ? 2 + 64 X ? 2 - 256L^ lX 2 2X3 

- 128i^ lX ? 2 +A(x P )n 16 [-1/24 XgXl2 - 1/8 Xl2X3 + 1/4 i?^ X i2 X3 + 1/12 ^ X L] 

- A( Xp )i5 [2 XpXl2 + 8 J RP XpXl2 + 2^ Xl2 ] - A( Xp )L r 3 [5 Xp Xi2 + 2R p q XpXi2 + 5R^xi2] 
+ 16 A( Xp ) LI RP X i2X3 + Mx P )Ll [4 XpXl2 + 8i^ XpXl2 ] - 16 A( Xp ) L r 6 RP X12X3 

+ A( Xp )L r 7 [16/3 R d Xp s + 32/3 i^ Xl2 ] ~ Mx P ) L r 8 [8 XpXl2 + 16/3i^ XlX2 + 32/3i^ X 2 2 + 4/3 (R d p ) 2 ] 

+ A( Xp ) 2 [5/32 Xl2 + 1/8 (i??) 2 X i2- 1/24 R*Xp] + Mx P )A(x p3 ) [l/24 Xp - 1/8 Xg + 1/24 i?£ Xp 

+ 19/72 RP Xq ] + A(x P )A( Xq 3) [1/4 X i 2 + 1/6 R p qXq ~ 1/24 R q pXp + 1/72 J2» X J 

+ A( Xp )A( Xl2 ) [1/18 Xp - 1/36 xia] - l/8A( Xp )A( X3 ) X3 + A( Xp )B( Xp , Xp , 0) [ll/24i?£ XpXl2 

- 1/12 ^ xl~ 1/8 fl?X?2 + 1/24 i^fc,] " A( Xp )B( Xq ,x q ,0) [l/24 RP q R d qXq + 1/16^x12] 
+ l/6^( Xp )S( Xl , X2 ,0)i?9 XpX i2 + l/12A( Xp )B( Xl , X2 ,0,fc)i?« Xp + l/4l( Xp! e)7r 16 [ XpX i 2 

- ^9X12X3] ~ ^(Xp3)7Tl6 [l/3 XpX 12 - 2/3 X i2 X3 - 1/3 X ? 2 ] + 8A( Xp3 ) J Lo Xp3X i2 
+ 20^( Xp3 )L3 Xp3X i2 - 16 A( Xp3 ) L£ Xp3X i 2 + 32^( Xp3 )Lg Xp 3 X i2 

+ l/4A( Xp3 )5( Xg , Xg ,0)^ Xl2 - A( Xp3 )S( X i, X2l 0) [4/9 Xp3X i2 + l/9 X i X2 + l/9i?« X 2] 

+ l( Xp3 )B(xi,X2,0,fe) [-2/9 Xp3 -2/9 Xl2 + l/9i?« X ,] - A( Xp3 ,e) 7T16 X12X3 

+ A(xi)A( X2 ) [-I/I6X12- 1/24x3 + l/4i^i2?xia] + 8A( Xl2 )I^ X 2 2 + 20 A( Xl2 ) U 2 x\ 2 

+ 32A( X i2)L r 6 xl 2 - l/4A( Xl2 ) 2 X i2 + A( Xl2 )S( Xl)X2 ,0) [1/9 XiX2 + 4/9 X \ 2 ] 

+ 4/9l( X i 2 )5( XllX2 ,0,fe) X i 2 - 2yl( Xl2 ,£)7r 16X 2 2 - A( Xl3 )l( X23 ) Xl2 + 48 A( Xs ) L\ X12X3 

+ 12A( X 3)L r 2 xi2X3 ~ 48A( X3 )L^ X i 2X3 + 48 A( X3 ) L r G X12X3 + 1/2 A(x 3 )B(xi, X2, 0) X12X3 

+ l/4A(x3)S(xi,X2,0,fe) X3 - 3/4A( X 3,e)7r 16X i2 X 3 - B( Xp , X P , «) Jr w [1/24^X9X12 + 1/8^X12X3] 

- 2S( Xp , Xp ,0)^^ XpXl2 - 5B( Xp , Xp ,0)L3i? pXpX i2 + 8B( Xp , Xp ,0)L4i?^ XpX i2 X3 

+ S( XpiXp ,0)Lg [4^ xp + 2^ X i X2 ] - 165( XpiXp ,0)^^ XpXl2X3 - B( Xp , Xp ,0)££ [8R p qX 3 p 

+ 4^ XlX2 ] + B( Xp , Xp ,0) 2 [9/32 RPR^XvXii- 1/32 R p q R d Xq Xi2] 

+ 1/ 6 B(x P , Xp , 0)B( X i, X2 , 0) R%Rp XpXi2 + l/12B( Xp , Xp ,0)B( Xl , X2 ,0,fc)i?^ Xp 

+ l/4B( Xp , Xp ,0,e)7r 16 ^ Xp3X i2 - l/16B(xi,Xi,0)B(X2,X2,0) R1R2X12 

+ 32/3B(xi,X2,0)L r 7 R d 1 R d 2X i2 + 16/3 B( X i, X2 , 0) L r s R d R d Xl2 + 16/3 S( X i, X2 , 0, /c) Ly i?^i? 2 
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+ 8/3B( X i,X2,0,k)L r s R d 1 R d 2 + H F (1, Xp , X P , X12, X12) [-1/4 Xp Xi2 - 15/32 X 2 2 + 3/8 i^xi 2 

- 1/8 (R P q ) 2 X 2 12 ] + H F (l, Xp ,Xl3,X23,Xl2) [1/4X12 + 3/4i^x P Xi2 + l/4i^X,Xi2-l/4i^Xi2X3] 

- F F (l,Xi,Xi2,X2,Xi2) [1/16 X12- 3/4 i?2^?X? 2 ] + l/4ff F (l,Xi2,Xi2,Xi2,Xi2)x?2 

+ 3/4tf F (l,xi 3 ,X23,X3,Xi 2 )xi2X3 - ff F (2, Xp ,X P ,Xi2,Xi 2 ) [l/8i2£ij£x? 2 -5/16i#x? 2 ] 

- !? F (2,x P ,Xi2,X 9 ,Xi2) [l/8i?^ x f 2 -l/16^ X y - l/4H F (2, Xp , Xl3 ,X23,Xi2)i?pX P 3Xi2 
+ 5/32i/ F (5,x P ,X P ,Xi2,Xi2)(i?p) 2 X? 2 + l/16i/ F (5,xi,X2,Xi2,Xi2)i?^ 2 X? 2 

+ ^f(l,X P ,X P ,Xi2,Xi2) [5/2xl 2 -R p q R q p X 2 i2] - 2#f(l,X P ,Xl3,X23,Xl2)i?2X?2 

- 2 J ff 1 F (l,X P 3,X 9 3,Xp,Xl2)X? 2 + H( (1,X12,X1,X2,X12) [1/2 X12 - ^2^1 X? 2 ] 

- l/2H[(3,xi2,X P ,Xp,Xi2) R q p R d p xi 2 + l/2H[{3,Xi2,Xp,X q ,Xi2)R q p Rtxi2 

- l/4i/ F (7,xi 2 ,X P ,X P ,Xi2)(^) 2 X? 2 - 3/4i/ F (l,x P ,X P ,Xi2,Xi2)x 2 2 
+ 3/4i? F i(l,x P ,Xi3,X23,Xi2)-RpXi 2 + 3/4^(1, 

X P 3 ? Xg3 ? X P 7 Xi2)i??X 2 2 

- 3/4ff F (l,x P 3,X 9 3,X 9 ,Xi2)i? P X 2 2 + ff2i(l,Xi2,X P ,X P ,Xi 2 ) [15 /32 X 2 2 - 3/ 8 R p R p xi 2 ] 

- i/ 2 F 1 (l,Xi2,Xi,X2,Xi2) [3/16 X 2 i2- 3/4 i? 2 i? 2 X? 2 ] + 3/4 if F (l, xi 2 , X12, X12, X12) X? 2 
+ 9/4 #£(1, X3,Xi3,X23,Xi2)x? 2 - 3/4^(3, x P 3, X P , Xg3, X12) -Rp X12 

- Hfi(3,xi2,X P ,X P ,Xi2) [3/8i?^x 2 2 ~3/16i?^x 2 2 ] - #£(3, X12, Xp, X„ X12) [3/8i?^x? 2 

- 3/16 R d pX i 2 ] + 3/32H F (7,xi2,X P ,X P ,Xi2)(i?^) 2 X 2 2 + 3/16 F (7, X12, Xi, X2, X12) i^xL (32) 

Cops = ^16^0 [8X7rXi 2 + 6x 2 Xi2 + 4xiXi2X2-4xi2X3X4-8x? 2 ] + A ^w L iXi2 +ni6L 2 [6x 2 Xi2 + 2x? 2 ] 
+ n le L r 3 [6x7rXi 2 +4x 2 Xi2 + 3xiXi2X2-3xi2X3X4-8x? 2 ] + "le [35/144 x^Xi 2 + 2 /3 x 2 Xi2 
+ 17/96 x1X12X2 - 17/96x12X3X4 - 17/72 x? 2 ] + 256L^ X 2 2 X34 - 512^x12X34 

- 256L^x 2 2 X34 + 128L^ lX 2 2 X34 + 256Lf X12X34 - 256L^x 2 2 X34 - 128 L r 5 L r 8 X l 2 

+ 64L5L; iX ? 2 + 64Lf X ? 2 - 256L^ lX 2 2 X34 - 128 L r 8 L r n X 3 12 + A( Xp )n 16 [l / 4 R p qn X12XU 
+ 1/12 RP q „ X \ 2 - 1/24 R% Xq Xi2 - 1/8 R c p X12X34] - A( Xp ) L£ [8 X P Xi2 + 2 R c p Xp Xi2 + 2 R d p xi 2 ] 

- A{ Xp )Ll [2^ 7r x P Xi2 + 5^x P Xi2 + 5^xi2] + 16A( Xp )L^xi2X34 + ^(x P )^5 [8^X P Xi2 
+ 4i?^x P Xi 2 ] - 16A( Xp )L^Xi2X34 +A( Xp )L r 7 [32/3 R d Xp + 16/3 R d Xq - 8/3 {R d f\ 

- A( Xp ) LI [16/3 iJ^ X1X2 + 32/3 R p qn X \ 2 + 8 R c p X P Xi2 + 4/3 (i?^) 2 ] 

+ A( Xp ) 2 [1/8 X12 + 1/8 2 Xi2 + 1/24 R p q7r R c p Xp - 1/24 (R c p ) 2 Xp + 1/32 (R c p ) 2 Xl2 ] 
+ A( Xp )A( Xps ) [1/48 i^Xp + 19/144 R p q7rXq ~ 5/144 RP S7rXp - 1/16 i^X, + 1/18 ^Xp] 
+ A( Xp )^(x 9S ) [1/12 Xq ~ 1/48 i?L^ sp X P + 1/144 R p S7r R z qsp Xq + 1/8 R c p X12] 

+ A(x P )A(x.) [1/24 i?^i?J p x P - 1/40 R^R^ Xq + 13/120 i?^i?? 2 Xp + 1/8 R p q7I R* 12 Xq ~ 1/60 R p ^ qq X « 

- 1/24 R* qq R c p Xp - 1/24 R c p Xn + 1/96 i$i£ 12 Xp + 1/32 R c p Rl 12 Xq - 1/16 R c p R v nl2 Xn - 1/12 R* 34p ] 
+ A( Xp )A( Xl2 ) [1/18 Xp- 1/36x12] - A(x P )A( X 3i) [l/12R p 7vX3i - 1/9 R; pXp ] 

+ A( Xp )B( Xpi Xp , 0) [H/24 R^R c p Xp Xi2 - 1/12 R p qiT R c p X % - 1/8 R p qn R c p x\ 2 + 1/24 R p qn R d x P ] 

+ A( Xp )B( Xp ,Xn,0) [l/6i?^ Xp + 1/12^2^X1X2] ~ A(x P )B(x qiXq ,0) [l/2ARP qn R d Xq 

+ 1/16 R c p R d X12] + l/6A(x P ) J B(xi,X2,0)i?^XpXi2 + 1/12 A( Xp )B(xi, X2, 0, fc) x P 

+ 1/4 A(x P ,e) Trig [XpXi2 - -R^ X12X34 - X12X34] + A(x P s)^ie [1/6 X9X12 + 1/4 X12X34 + 1/12 Xi2Xs 

- 1/6 Xi 2 ] + 4:A(x P s)LqX P sXi2 + 10A(x P s)L 3 x P sXi2 ~ & A(x P s) L r 5 x P sXi2 + 16 A(x ps ) L r & x P sXi2 
+ A(x ps )A{x.) [1/18 R; p x p - 5/144 i?^x P - 1/16 i^Xg + l/8i2J,Xi2 - 5/48^ sXp 

- 19/144 R^ qs R z ps7r Xq + 1/48 R\ 2 Xp ] - A( Xps )B(x P , Xn,0) [2/9 R% t Xp X P s + 1/9 Rl s x P Xi2 
+ l/9i^ s x 9 X ps + l/18i?^ ST xiX2 + l/18i?r 2 X P ] + 1/8 A(x ps )B(x q ,X q , 0)^X12 

- A( Xps )B(xi,X2,0) [1/18 i?^ Xp + 2/9 i?Lx P ,Xi2 + 1/18 i?LxiX 2 ] 
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+ A{ Xps )B( XuX2 ,0,k) [1/18 R% Xq - 1/9 R q S7T X P s- 1/9 R q S7r Xi2] - 1/4 A{ Xps , e) tt 16 [ Xl2X3 4 + XuXs] 

- l/AA{x P z)A{ Xq i)xi2 - l/12A(x,r)7ri6 [K12 + K12 X12] 

- A( Xw )LS [12^2X^X12 + 2^12X^X12] + 16A(x,)LIXirXi2 + 4 A( X .) £2X^X12 

- A( Xir )LS [12^2X^X12 + 5^12X^X12] - 16A( Xw )L5 [XttXi2 -^2X^X12] 
+ A( X7r )L^ [4i?; p x w Xi2+8i?? 2 x? 2 ] + 16A( X7r )Lg [XttXi2- ^2X^X12] 

+ A( Xn ) L r 7 [64/3 i?? 2 X.X12 - 64/3 R\ 2 X \ 2 + 16/3 ^ 3 4 P X q - 8/3 (^4)^12] 

- A( x .) LI [32/3 i?? 2 X1X2 + 64/3 i?? 2 x? 2 + 16/3 < 12 x 2 + 4 < 12 X1X2 - 4/3 R v nl2 X3Xi ] 
+ A( Xjr ) 2 [1/8 (i?? 2 ) 2 xi2 - 1/24 Rl 2 Rl 12 x. + 1/24 i^ 2 Ki2 X i 2 - 1/24 R v nl2 Xn 

- 1/24 (Ki 2 ) 2 X. + 1/32 « 12 ) 2 Xi 2 ] - A(x*)A(x34) [1/3 R n 12 X12 + 1/36 R v nl2 Xn ] 

+ A(x„)B(x P ,Xp,0) [l/4R p qv RZ P X 2 p - 1/12 R; p R d pXp - 1/24 R* qq R d Xp + 1/96 R d R v nl2Xp 

+ 1/32 R d p Rl l2 Xq - 1/16 R d R v nl2 X,} + M.X*)B(Xp, Xn,0) [l/6 i^i?? 2 x P Xi2 + 1/12 i?^ 2 Xp] 

- l/4A(x.)B(x.,X.,0)i?^ 2 xi2X34 + 1/6 A( X7r )B(xi,X2, 0)^^X12X34 

+ l/12A{ X ^B{ Xl , X2 ,0,k)R^R 2 ^ X34 + A{ X ^,e)n 16 [-1/4 X .Xi2 + 1/4^2 X-X12 + 1/4 K12X.X12] 
+ A(xi)A( X2 ) [- 1/24 R p q7r R c qXq + 1/4 R^R^ Xl2 - 1/16 ^ X u] + 8i(xi 2 )^X 2 2 
+ 20A(xi 2 )^X 2 2 + 32A(xi 2 )^X 2 2 - l/4i(xi 2 ) 2 Xi2 + i(xi 2 )B(xi, X2, 0) [l/9 X1X2 + 4/9 X ? 2 ] 
+ 4/9A(xi 2 )B(xi,X2,0,fc)xi2 - 2A(xi2,£)7Ti 6 x 2 2 + 32 i( X34 ) ^ X12X34 + 8 A(x 34 ) L r 2 X i2 X 34 

- 32A(x 34 )L4Xi2X34 + 32A(x 34 )^Xi2X34 - A( X34 ) B ( Xp , X71 , 0) [2/9R^ 2XpX7I + l/9Rl 2Xq x, 

+ 1/6^ 2 X1X2] + l/2A( X 34)S(X7r,X7r,0)i?r 2 XxXl2 + A( X34 )B( X 1 , X 2 , 0) [1/3X12X34 + 1/9^X1X2 

+ 1/3^2X12X34- 2/9 Rl 2X \ 2 ] + A( X34 )B( XuX2 ,0,k) [1/6 X34 + 1/9 i?? 2 X i 2 + 1/6 i?? 2 X34 ] 

- 1/2^4(X34,£)7T16X12X34 - l/4A(xi 5 )A(x 2s )xi2 - B( Xp ,Xp,0)tt 16 [l/24R d XgXl2 + l/8i?pXl2X34] 

- 2B(x P ,Xp,0)L r R d X pXi2 - 5B( Xp , Xp ,0)L r 3 R d XpXl2 + 8B( Xp , Xp ,0)L r 4 R P 7rXpXl2X34 
+ B( Xp ,Xp,0)L r 5 [4i?^ X p + 2^X1X2] ~ 16S(x„,X P ,0)^i?^x P Xi2X34 

- B( Xp , Xp ,0)L r 8 [8R p qwX 3 p + 4^XiX2] + B(xp,X P ,0) 2 [9/32 R? q „R d Xp Xi2 - 1/32 RP qn R d XqXl2 ] 

+ S(xp,Xp,0)B(x P ,X^0) [l/6K 2 i?^x 2 + l/12i?r 2 i?^XiX 2 ] + 1/6 S(x P , Xp, 0)S( X i, X2, 0) XpX i 2 
+ l/12B(x P ,Xp,0)B(xi,X2,0,fc)i?^x P + l/8B(xp,Xp,0,e)7r 16 [i^ x P Xi2 + R$ X12X34] 
+ B(x p ,Xx,0)L; [16/3 R d R z 7r34qXp + 32/3 i^i?* 34g X12] + B(x P , X., 0) Lg [8/3^ 34gX p 
+ 16/3i?^ 34(? xi2] + 8B( X;r ,x.,0)^i?r 2 x 2 Xi2 + 4B( X;r ,x 7 r,0)^i?r 2 x 2 Xi2 

- 16S( X7r ,x 7r ,0)^^ 2 x 2 Xi2 - 16 5( XT;XT ,0)L^ 2 (i?* 4 ) 2 Xl2 - 4B( Xn , Xw ,0) L r 8 [Rl 2Xl2X \ 
+ R12 Xi2xl] - l/16B( Xl , Xl ,0)B( X2 , X2l 0)R d R d 2Xl2 + 32/3 B( Xl , X2 , 0) V 7 R d R d 2 Xl2 

+ 16/3B( X i,X2,0)^^ J R 2 d Xl2 + 16/3 B( Xl , X2 ,0, A) + 8/3 B( Xu X2 , 0, fc) R d R d 

+ H F {l, XpiXpiXl2lXl2 ) [-l/4 XpXl2 - 3/8 x? 2 - l/8(^J 2 x 2 2 + 3/8i?^x 2 2 - 3/32 (R c p ) 2 X \ 2 ] 

+ H F (l, Xp , xi„ X2s, X12) [-5/8 R p q7T XqXl2 + 1/4 RP qn XpsXl2 + 3/16 RP S7r XqXl2 + 7/16 R p „ Xl2Xs 

- 1/8 Rp XpsX 12] + l/4ff F (l,Xp3,X 9 4,X34, X 12)Xl2X34 + H F (1 , x ^ , Xp , Xl2 , X 12 ) [l/2 iS^iZJa X?2 

- 1/8 i^iCi2X 2 2 + 1/4 V^?2 + 1/16 i?pKi 2 X 2 2 ] + ff F (l,X.,X.,Xi2,Xi2) [l/2(i?? 2 ) 2 X 2 2 

+ 1/2 K 2 Ki 2 X 2 2 + 5/32 (Ki 2 ) 2 X 2 2 ] + ^ F (l,X^Xis,X2 S ,Xi2) [-3/16 R" ps XpXl2 - 5/16 i?; s X? Xi2 
+ 1/2 i?r 2 x 2 2 - 1/16 KpsXtXu- 1/16 KuXuXs] + H F (l, Xl , Xl2 , X2 , Xl2 ) [-1/8 R p qn R q X \ 2 
+ 3/4Rl w R 2 ln x 2 2 + 1/16 R\R C 2X \ 2 \ + l/4H F (l lXl2 , Xl2lXl2 , Xl2 ) x 2 12 

- H F (2 lXp , Xp , Xl2lXl2 ) [1/8 R p qn R d pX 2 2 - 5/16 R c p R d X 2 2 ] + H F (2, Xp , Xn , Xl2 , Xl2 ) [l / 4 R; p R d X 2 2 
+ l/16i?Xi2X 2 2 ] - i? F (2,Xp,Xi2,X 9 ,Xi2) [l/8i?^x 2 2 -l/16i?^pX 2 2 ] 

- l/8H F (2,Xp,Xis,X2 S ,Xi2)#pXp S Xi2 + 5/32i/ F (5,Xp,Xp,Xi2,Xi2)(i?p) 2 X 2 2 

+ l/16H F (5,xi,X2,Xi2,Xi2)i?fi? 2 d X 2 2 + #r(l,Xp,X P ,Xi2,Xi2) [2x 2 2 + KJ 2 xL-^^pX 2 2 
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+ V2 ORp) 2 Xi 2 ] + H[ (l, Xp ,Xis,X2s,Xi2) R p q7T R z sqp xl 2 ~ H[(l, X ps,Xqs,Xp,Xi2)RLxl2 

~ H 1 {1, Xps, Xqs, Xir, Xl2) RuRspirXu + 1/ 2 #1 ( 1 , \\2 , Xp ) Xtt , Xl2 ) [Rq^r Rl 12 Xl2 ^ RppRqp^Rp X12] 

- Hf(l,xia,X*,X*,Xi2) [l/2i?Mi2X?2 + l/4(Ki2) 2 X? 2 ] +^f(l,Xi2,Xi.Xa,Xia) [l/2i?^ X ? 2 

7 Xl2) Xpi Xpi X12) [-Kqir-KpXl2 RpRp X12] 

+ l/2fff , (3,Xi2 J X P ,X g ,Xi2)^ r i^Xia + 1/2 #f (3, X12, X P , X., X12) R^Ki^t ^ 

- l/4H[ (7,Xl2,Xp,Xp,Xl2) (Rp) 2 Xl2 - 3 / 4 - ff 2l( 1 'Xp,Xp ; Xl2,Xl2)Xl2 

- 3/8i/ 2 F 1 (l,X p ,Xl S ,X2 S ,Xl2)^ 7r -RsgpX?2 + 3/8if^(l,Xp S ,Xg S ,Xp,Xl2) [^Xl2+CXl2 

- RpXu] + 3 / 8 ^2l( 1 >Xp«)XgS)Xg ! Xl2) Rp^R Z S pqXl2 + 3/8 Hili 1 ' Xps, Xqs, Xtt, Xl2) R^RpqnRlpnXu 

~ 3/8H[ 1 (l,x„Xi s ,X2s,Xi2)Ri 2 R Z si,R Z s2,xl2 + fl£(l,Xi2,Xp.Xp.Xi2) [3/8 X ?2 + 3/8(i?^) 2 x? 2 

- 3/8 ^^Xi2 + 3/32 (i?p 2 x? 2 ] - Hf 1 (l,Xi2,%,,X ) r,Xi2) [ 3 /8 ^Ki2 X? 2 - 3/16 R c p R v nl2 Xu] 

+ 3/32^ 2 Jr 1 (l,Xl2,X7r I X 7 r,Xl2)(Kl2) 2 X?2 + ff 2 F l(l,Xl2,Xl,X2,Xl2) [-3/8 X 12 + 3/4 Xl2 

+ 3/16^^2X12] + 3/4if^(l,Xl2,Xl2,Xl2,Xl2)Xl2 + 3/4i/ 2 F 1 (l,X34,Xp3,Xg4,Xl2)Xl2 

- 3/81f 2 F 1 (3,x ps ,Xp,X g5I Xi2)^pX?2 - 3/8£f£(3,xi2,Xi»,Xi.,Xia) fei?pX 2 2 - l/2i$J#Xi a ] 

- Hfi(3,xi2,%,>X„Xi2) [3/8i?^x 2 2-3/16i?^ X 2 2 ] + 3/16 ff£(3, Xl2 , X P , X., X12) R$Ku X? 2 

+ 3/32i/ 2 F 1 (7,xi2,X P ,X P ,Xi2)(i?^) 2 X 2 2 + 3/16i7 2 p 1 (7,xi2,Xi,X2,Xi2)i?^X 2 2 - (33) 



V. NUMERICAL ANALYSIS AND 
CONCLUSIONS 

The NNLO results given in the preceding section de- 
pend on a number of 0(p 4 ) and 0(p 6 ) LEC:s. Eventually, 
their proper values, from which the physical LEC:s l\ and 
c\ can be obtained, should be determined by a fit of the 
NNLO expressions to lattice QCD data. For the present 
purpose, in order to have some reasonable estimates of 
the LEC:s, we simply use the values obtained from ex- 
periment in Ref . | l5| , which used three-flavor continuum 
XPT to NNLO. 

The combination of LEC:s used in this paper corre- 
sponds to 'fit 10' of Ref. pi, which had F = 87.7 MeV 
and fx — 770 MeV. The parameters which were not de- 
termined by that fit, namely Lq, L\ and Lg, have been 
set to zero at the abovementioned scale n, as have all the 
KJ. Also, L\ x has been treated according to the condi- 
tion L r u = -III A = -(2L r 4 + LI). It should be noted 
that Lq cannot be determined from experiment, but is 
obtainable from partially quenched simulations. Some 
recent results on L\ and Lg may be found in Ref. [l6j| . 
If we temporarily reintroduce the superscripts (2pq), the 
choice of LEC:s for the numerical analysis is 



L 



F 


= F 


fit 10 ' 




r(2pq) 
'i=1...10 


= L\= 


1...10 


fit 10 ' 


r r(2pq) 
^11 


= ~(2L r 4 + 


LI) 


fit 10 ' 


T r(2pq) 


T r{2pq) 
- "12 - 


j^r{2pq) _ q 



In order to set the scale for the lowest order masses, for 
a quark mass of 100 MeV we have roughly x ~ 0.01 GeV 2 
and for 500 MeV we have x ~ 0.25 GeV 2 . A crude 
estimate of the range of validity of PQxPT is for all 
XnXij ~ 0.3 GeV 2 . 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
Xi [GeV 2 ] 



FIG. 3: The NNLO modification Af of the decay constant for 
dvai = 1 and (tea = 1 as a function of the valence quark mass 
parameter xi an< l the sea-quark mass parameter \3- The 
function plotted represents the sum of the 0(p 4 ) and 0(p 6 ) 
contributions. 
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FIG. 4: Comparison of the NLO and NNLO corrections Af to the decay constant for d va i = 1 and d sea = 1. The curves labeled 
'NLO' represent the 0(p 4 ) contribution to the decay constant, and those labeled 'NNLO' represent the sum of the 0(p 4 ) and 
0(p b ) contributions. The dependence on the sea-quark mass X3 has been parameterized as \3 = Xitan(#). Present Lattice 
QCD simulations are performed in the region where 9 > 45° . 



A. Numerical Results for Decay Constants 

The modification of the decay constant, when correc- 
tions up to NNLO are taken into account, is parameter- 
ized in terms of the relative change A p , which is defined 

as 

F phys = F (1 + A F ) , (35) 

so that Ap is expected to disappear in the chiral limit, 
when both sea and valence quark masses are set to zero. 
That this condition is satisfied can be seen directly from 
Fig. [3J However, if the sea-quark mass is set to zero 
for a constant valence quark mass, the decay constant 
approaches a constant value. This behavior is expected 



and may be attributed to the unphysical quenched chiral 
logarithms and similar effects. 

The numerical results for the pseudoscalar meson de- 
cay constants are shown in Figs. |3] and 21 For simplicity, 
we have considered the case with d va i = 1 and d sca = 1, 
i.e. the valence quark masses and sea-quark masses are 
degenerate, but different from each other. Thus the re- 
sults presented correspond to Eq. (|16fl for NLO, and 
Eqs. (O and (JTHJ) for NNLO. 

The convergence of the chiral expansion for the decay 
constant is investigated in Fig. 0] for different combina- 
tions of sea and valence quark masses. We have chosen 
the sea-quark mass larger than the valence quark mass 
in our examples since this is the most likely regime for 
lattice QCD simulations. Overall, Fig. [3] shows that the 
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FIG. 5: Comparison of the NLO and NNLO corrections Am to the meson mass for d va i = 1 and d sea = 1. The curves labeled 
'NLO' represent the C(p 4 ) contribution to the meson mass, and those labeled 'NNLO' represent the sum of the 0(p 4 ) and 
C(p 6 ) contributions. The dependence on the sea-quark mass X3 nas been parameterized as X3 = Xitan(#). Present Lattice 
QCD simulations are performed in the region where 9 > 45°. 



NNLO corrections to the decay constant are significant 
at least for our choice of LEC:s. Nevertheless, the com- 
parison between the NLO and NNLO curves in Fig. 0] 
also shows that there exists a rather large range of in- 
put quark masses for which the convergence of the chiral 
expansion is reasonable. The cases with less degenara- 
cies in the input quark masses could also be considered, 
but such plots are difficult to present in an instructive 
way since these functions depend on up to four different 
quark masses. One possible way to proceed is to look at 
various special cases where the ratio between the sea or 
valence quark masses is fixed. To illustrate the effects 
of nondegenerate sea and valence quarks, such plots are 
given in the next subsection for the meson mass. 



B. Numerical Results for Meson Masses 

The numerical results for the pseudoscalar meson 
masses are shown in Figs. [5] through [S] Part of the pre- 
sentation is similar to that of the decay constant results 
in the preceding section, i.e. results are given for c? va i = 1 
and d sca = 1. These results correspond to Eq. (|2"fj)l for 
NLO, and Eqs. iJUJ and (gSJ for NNLO. In addition some 
results are also presented for c? va i = 1, c? sca = 2 in Fig.0 
and for d va \ = 2, tf soa = 1 in Fig. El The modification of 
the meson mass when NNLO corrections are taken into 
account is parameterized in terms of the quantity Am, 
which is defined as 



Mph ys = Ml (1 + A M ) , 



(36) 
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so that Am represents a relative deviation from the low- 
est order result Mg . With this definition, Am is expected 
to disappear in the chiral limit, when both sea and va- 
lence quark masses approach zero at a constant ratio. 
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-0.1 
-0.2 
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FIG. 6: The NNLO modification Am of the meson mass for 
dvai = 1 and d S ea = 1 as a function of the valence quark mass 
parameter Xi an d the sea-quark mass parameter X3- The 
function plotted represents the sum of the 0(p 4 ) and 0(p 6 ) 
contributions. Note the divergence of the result near xi = 0, 
which vanishes for X3 = 0. 



In general, the appearance of Fig. [S] indicates that 
the NNLO contribution to the meson mass is more pro- 
nounced than for the decay constant, which is due to 
the fact that the NLO contribution to the meson mass 
vanishes for certain values of the input quark masses. 
Another noteworthy feature of Fig. is the logarithmic 
divergence for small values of the valence quark mass xi- 
This feature, which is well-known at NLO and persists 
at NNLO, is due to the unphysical quenched chiral loga- 
rithms. This divergence disappears if the sea and valence 
quark masses are sent to zero in constant ratio, which is 
also confirmed by Figs. [S] through |SJ 

There is no a priori reason to believe that the con- 
vergence of the meson mass should be intrinsically worse 
than that of the decay constant. Nevertheless, the com- 
bined NLO + NNLO result in Fig. [B] vanishes in several 
places, with the NLO and NNLO contributions canceling 
each other. The decay constants showed a more uniform 
convergence for the masses and input LEC:s used here. 
In this respect the curves shown in Fig.[S]are less instruc- 
tive as to the convergence of the chiral expansion than 
those shown in Fig. 0] Admittedly, however, the inclu- 
sion of NNLO effects into the pseudoscalar meson mass 
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FIG. 7: The modification Am of the NNLO meson mass for 
dvai = 1 and d sca = 2. The sea quark mass dependence has 
been parameterized as Xi ~ %X'a an d the relation between the 
sea and valence quark masses as xz — Xi tan(#). 
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FIG. 8: The modification Am of the NNLO meson mass for 
dvai = 2 and d SC a = 1. The valence quark mass dependence 
has been parameterized as X2 = xxi and the relation between 
the sea and valence quark masses as X3 — Xi tan(#). 



in PQxPT is a major effect which changes features of the 
NLO calculation. 

It is also noteworthy that the convergence of the chiral 
expansion for the quantities considered apparently de- 
pends rather strongly on the values of the LEC:s. For 
example, the obvious trial choice of setting all the L\ to 
zero leads, incidentally, to a much larger NNLO contri- 



19 



bution which converges less well also for the decay con- 
stants. Setting several of the K? to nonzero values can 
also change this picture quite dramatically. Thus a defi- 
nite statement about convergence can only be made when 
some estimates are available for the K\. 

C. Determination of LEC:s from Lattice QCD 
Simulations 

In the previous sections, a sample of numerical results 
have been presented in order to illustrate the effects of the 
NNLO contributions when the LEC:s are known. Since 
the intended area of application of the NNLO formu- 
las given in this paper is the determination of LEC:s 
from PQQCD Lattice simulations, some attention should 
also be devoted to finding efficient methods for this task. 
Again, it is understood that the jj^ 1 ^ and K^ 2pq ^ are 
simply denoted by L\ and K\. 

At NLO in PQxPT, the expressions for the pseu- 
doscalar meson masses and decay constants can be 
cleanly separated into a tree-level contribution which is 
a function of the L\ and a remaining part, independent 
of the L\, which involves the chiral logarithms. Thus at 
NLO it is possible to develop various optimal strategies 
for the purpose of obtaining estimates of the L\ by fits 
to Lattice QCD data £|. At NLO, the dependence on 
the L\ is proportional to 2L\ X34 + £5 X12 for the de- 
cay constant and (4Ljj — X34 + (2Xg — L£) X12 for 
the masses, from which an overall factor X12 has been 
removed. The LEC:s can thus be determined from the 
1+1 case by using a constant sea quark mass and vary- 
ing the valence quark mass. The more complicated cases 
with nondegenerate quarks, i.e. xi 7^ X2 or X3 + 1 X4 are 
not needed at NLO. 

A similar analysis which involves the determination of 
LEC:s from a fit to Lattice QCD simulations using NNLO 
PQxPT is more challenging, not only because of the 
complexity of the NNLO expressions, but mainly since 
the chiral logarithms now involve the L\. Furthermore 
the coefficient functions of the L\ at NNLO have a more 
complicated dependence on both sea and valence quark 
masses. On the other hand, the situation for the K\ at 
NNLO is similar to that of the L\ at NLO, although the 
much larger number of K\ parameters, along with the 
appearance of contributions bilinear in the L\ certainly 
complicates the analysis. 

In general, at 0(p 6 ) the expressions for the pseu- 
doscalar meson masses and decay constants depend on 
12 K\ parameters, 5 for the decay constants and an ad- 
ditional 7 for the masses. The dependence of the decay 
constant on the K\ is proportional to 

16X 12 ^[ 9 " 8x1X2 (K r 19 -K r 23 ) + 16xi2X34^ 2 r o 
- 16 X3X4 K r 21 + 32 xIa (K r 2 i + K r 22 ) , (37) 

from which all 5 appearing constants, K[ 9 23 can be de- 
termined but at least one case with nondegenerate sea 



quark masses needs to be calculated. From the 1+1 case, 
three combinations can be determined. Also, to separate 
K\ 9 from K 23 a case with nondegenerate valence quark 
masses, d V ai = 2, is needed, otherwise only the combina- 
tion K[q + K 23 can be obtained. Likewise nondegenerate 
sea quark masses, d soa = 2, are needed to disentangle K 21 
from K 22 , simulations with degenerate sea quark masses 
allow only the determination of K 21 + 2K 22 . Similarly, 
the dependence of the meson masses on the K\ is 

- X12 (32 K[ 7 + 32 K[ 9 - 96 K r 25 - 32 K r 39 ) 

+ X1X2 (16XT 9 -16*23 - 48KJ B ) 

- X12X34 (64 K{ & + 32 K r 20 - 64 K r 26 - 64 K r 4Q ) 
+ X3X4 (32^-32^) 

- X34 (64 K r 21 + 64 K r 22 - 64 K r 26 - 192 K r 27 ) , (38) 

from which an overall factor of X12 has again been re- 
moved. From the 1+1 case we can obtain 3 combinations, 
one more from a case with <i va i = 2 and one more from a 
case with d sea = 2. The particular combinations can be 
easily read off from Eq. 1)38(1. Thus, in order to obtain 
independent information on as many of the K^-.s as pos- 
sible, simulations with nondegenerate quark masses are 
clearly preferrable. On the other hand, a possible strat- 
egy for the determination of the L\ from NNLO fits is to 
use only one distinct sea and valence quark mass, since 
then the number of distinct combinations of K\ parame- 
ters is minimal, while the possibility of distinguishing all 
of the L\ 6 and Lg is retained. 

At NNLO, additional Lt are present in the expressions 
for the decay constant and the meson mass, namely L 7 
and the Lg 3. Of these, L 7 can be determined from the 
diagonal mesons Q by consideration of the coefficient of 
the double pole. This calculation is not yet performed 
up to NNLO. The other four L\ are those relevant at 
NLO for meson-meson scattering. Their values have to 
be used as input but could in principle be separated from 
the K\ since they multiply nonanalytic dependences on 
the quark masses. A detailed study of this is beyond the 
scope of the present work. 

The present analysis shows that while the number of 
free parameters in a fit of the NNLO expressions to Lat- 
tice QCD data is large, the possibility of producing a 
large amount of simulation data with different combina- 
tions of sea and valence quark masses should compensate 
for this. Furthermore, it should be kept in mind that 
many of the in principle free parameters L\ are already 
quite accurately known from previous work. 

D. Conclusions 

In summary, we have calculated the pseudoscalar me- 
son masses and decay constants to NNLO in Partially 
Quenched Chiral Perturbation Theory for two flavors of 
sea-quarks. Explicit analytical formulas of reasonable 
length and complexity have been given, along with a nu- 
merical analysis which demonstrates the numerical im- 
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plementation of the results. Although the NNLO results 
appear to be of a reasonable magnitude, definite state- 
ments concerning the convergence of the chiral expansion 
have to be postponed, as the magnitude of the results de- 
pends rather sensitively on a number of largely unknown 
LEC:s. 

An implementation of the present results will be made 
available by the authors in the near future 17]. The 
availability of these results makes it possible to perform 
a NNLO chiral extrapolation from partially quenched lat- 
tice data. However, usage of the present results requires 
Lattice QCD data on the pseudoscalar meson masses and 
decay constants which has been extrapolated to the con- 
tinuum limit and to infinite volume at different values 
of the quark masses. Although such data are not yet 
available, this situation will undoubtedly change in the 
near future. The ultimate goal of such fits to lattice data 



is thus to provide reliable estimates of the (unphysical) 
low-energy constants of PQxPT, from which the values 
of the physical ones can then be extracted. 
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